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Gas giant planets play a fundamental role in shaping the orbital architecture of planetary systems and 
in affecting the delivery of volatile materials to terrestrial planets in the habitable zones. Current theories 
of gas giant planet formation rely on either of two mechanisms: the core accretion model and the disk 
instability model. In this chapter, we describe the essential principles upon which these models are built 
and discuss the successes and limitations of each model in explaining observational data of giant planets 
orbiting the Sun and other stars. 



1. INTRODUCTION 

Jupiter and Saturn are composed predominantly of Hy- 
drogen and Helium and are therefore referred to as gas gi- 
ants, although most of these elements are not in gaseous 
form at the high pressures of these planets' interiors. The 
majority of extrasolar planets (or exoplanets) discovered so 
far have masses in excess of about one quarter of Jupiter's 
mass (Mj) and are known to be or suspected of being gas 
giants. Since Helium and molecular Hydrogen do not con- 
dense under conditions typically found in star forming re- 
gions and in protoplanetary disks, giant planets must have 
accumulated them as gasses. Therefore, giant planets must 
form prior to the dissipation of protoplanetary disks. Opti- 
cally thick dust disks typically survive for only a few mil- 
lion years and protoplanetary disks lose all of their gaseous 
contents by an age of < 10 7 years (see the chapter by 
Roberge and Kamp), hence giant planets must form on this 
timescale or less. 

Giant planets contain most of the mass and angular mo- 
mentum of our planetary system, and thus they must have 
played a dominant role in influencing the orbital properties 
of smaller planets. Gas giants may also affect the timing and 
efficiency of the delivery of volatile materials {Chambers 
and Wether ill, 2001), such as water and carbon compounds, 
to the habitable zones of planetary systems, where liquid 
water can exist on a (rocky) planet's surface. In the solar 
system, Jupiter is also believed to have reduced the impact 
rate of minor bodies (such as comets) on the Earth (Horner 
et al., 2010). Therefore, understanding giant planet forma- 
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tion is essential for the formulation of theories describing 
the origins and evolution of terrestrial planets capable of 
sustaining life in the form of complex organisms. 

Observations of exoplanets have vastly expanded our 
database by increasing the number of known planets by 
over one and a half orders of magnitude. The distribution 
of observed exoplanets is highly biased towards those ob- 
jects that are most easily detectable using the Doppler ra- 
dial velocity and transit photometry techniques (see chap- 
ters by Lovis and Fischer and by Winn), which have been 
by far the most effective methods of discovering planetary- 
mass object orbiting other stars. Although these extraso- 
lar planetary systems are generally different from the so- 
lar system, it is not yet known whether our planetary sys- 
tem is the norm, quite atypical, or somewhere in between. 
Nonetheless, some unbiased statistical information can be 
distilled from available exoplanet data (see the chapter by 
Cumming). 

The mass distribution function of young compact ob- 
jects in star forming regions (e.g., Zapatero Osorio et al., 
2000) extends down through the brown dwarf mass range 
to below the deuterium burning limit (12-14 Mj). This ob- 
served continuity implies that most isolated brown dwarfs 
and isolated high planetary-mass objects form via the same 
collapse process as do stars. 

However, star-like direct quasispherical collapse is not 
considered a viable mechanism for the formation of Jupiter- 
mass planets, because of both theoretical arguments and ob- 
servational evidence. The brown dwarf desert, a profound 
dip over the range from ~ 15 Mj to ~ 60 Mj in the mass 
distribution function of companions orbiting within several 
AU of Sun-like stars, suggests that the vast majority of gas 
giants form via a mechanism different from that of stars. 
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A theory based on a unified formation scenario for rocky 
planets, ice giants much as Uranus and Neptune, and gas gi- 
ants is the core accretion model, in which the initial phases 
of a giant planet's growth resemble those of a terrestrial 
planet's. The only alternative formation scenario receiving 
significant attention is the disk instability model, in which a 
giant planet forms directly from the contraction of a clump 
of gas produced via a gravitational instability in the proto- 
planetary disk. 

In this chapter, we introduce the basic physical concepts 
of gas giant planet formation according to core accretion 
models (Section 2) and disk instability models (Section 3). 
In Section 4, we present theoretical arguments suggesting 
that most of the giant planets known to date formed via core 
accretion and pose some of the outstanding questions that 
still need an answer. Future prospects that may settle some 
of the lingering issues are discussed in Section 5 . 

2. CORE ACCRETION MODELS 

The initial phases of a giant planet's growth by core nu- 
cleated accretion proceed through an accumulation process 
of solid material, in the same fashion as terrestrial planets 
form. Dust and small solid grains, which are entrained in 
the predominant gas component of a protoplanetary disk, 
coagulate into larger particles. Centimeter- sized particles 
tend to settle towards the disk midplane, aggregating and 
eventually forming kilometer- sized agglomerates, referred 
to as planetesimals. 

Planetesimals grow larger via pairwise collisions, lead- 
ing to the formation of a planetary embryo. An embryo may 
have the ability to grow at a rate that increases as its mass 
increases, eventually consuming nearly all planetesimals in 
the neighborhood, while rapidly gaining mass and becom- 
ing a planetary core or protoplanet. Once a planetary core 
grows large enough and the escape velocity from its surface 
exceeds the thermal speed of the surrounding gas, a tenuous 
envelope of gas begins to accumulate around the solid core. 

For most of the following planet's growth history, ther- 
mal pressure effects within the envelope regulate the ac- 
cretion of gas. The ability of the envelope to radiate away 
the gravitational energy released by incoming planetesimals 
and by contraction limits the amount of gas that can be ac- 
creted by the planet. A slow contraction phase ensues in 
which the accretion rates of both solid material and gas are 
small. However, as the protoplanet grows and its total mass 
exceeds the value beyond which the pressure gradient (in 
the envelope) can no longer balance the gravitational force, 
the envelope undergoes a phase of rapid contraction, which 
allows more gas to be accreted. The augmented mass of the 
envelope triggers further contraction and gas can be thereby 
accreted at an ever increasing rate. 

During this epoch, known as runaway gas accretion 
phase, the gas accretion rate is regulated by the ability of 
the surrounding disk to supply gas to the planet's vicinity. 
In this stage, once the planet's mass exceeds a few tenths of 
Jupiter's mass (for typical temperature and viscosity condi- 



tions in a protoplanetary disk around a solar-mass star), the 
gas accretion rate decreases as the planet mass increases due 
to tidal interactions between the planet and the disk. The gi- 
ant planet continues to grow at an ever decreasing rate until 
there is no gas available within the planet's gravitational 
reach. 

2.1. From Dust to Planetesimals 

The formation of a heavy element core is an essential 
part in any core accretion model. Hence, for completeness, 
we present a summary of the basic elements involved in 
this process. An in-depth discussion of these concepts can 
be found in the chapter by Chambers. 

The core formation starts at some distance from the star, 
most likely beyond the "snow line" where disk tempera- 
tures allow for condensation of water ice and solid material. 
This process begins from (sub-)micron-sized dust particles, 
which may have originated from the interstellar medium 
and/or condensed within the disk's gas. Such small solid 
particles are well-coupled to the gas, on account of their 
large surface area to mass ratio, and therefore move with it. 

As particles collide and stick together, they can grow 
larger and start to de-couple from the gas and interact with 
it because of differential rotation between solids and gas. 
The relative velocity stems from the fact that the gas ro- 
tates about the star slower than do solids due to the radial 
pressure gradient that partially counteracts the gravitational 
attraction of the gas towards the star. The resulting inter- 
action can be described in terms of friction via gas drag on 
the particle, producing an acceleration proportional to and 
in the opposite direction of the relative velocity between the 
particle and the gas. This friction generally removes orbital 
angular momentum from solid particles, causing them to 
drift towards the star. 

Along with the radial drift, small solids also experience a 
friction in the direction perpendicular to the disk midplane. 
Assuming that the vertical motion of the gas is negligible 
(compared to that of solids), the vertical motion of a solid 
particle is obtained by applying the second law of dynamics 

d 2 z 1 dz ~ 

in which z represents the distance over the disk midplane, 
Tf is the friction timescale (also referred to as "stopping 
time"), and ft is the particle's angular velocity about the 
star. The second term on the lefthand side of equation (1) 
is the opposite of the frictional force per unit mass exerted 
by the gas and the third term is the opposite of the vertical 
gravitational acceleration imposed by the star. In the limit 
Tf — y 0, dust is perfectly coupled to gas and dz/dt is equal 
to the vertical gas velocity, which is zero by assumption. 
Hence z is a constant and no sedimentation would occur. 
Notice that for the gaseous part of the disk, equation (1) is 
replaced by an equation for hydrostatic equilibrium in the 
vertical direction, with the vertical component of the gas 
pressure gradient balancing the vertical gravitational force 
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(see equation (3)). In the limit r/ — » oo, equation (1) as- 
sumes the form of a harmonic oscillator, hence particles os- 
cillate about the disk midplane with a period equal to 2tt/Q . 

For solid particles whose size is shorter than the mean 
free path of gas molecules and whose velocity relative to 
the gas is slower than the gas sound speed, c, the friction 
time is (Epstein, 1924) 



9 



(2) 



where p sp and R sp are the solid particle's volume density 
and radius, respectively, while p is the gas density. In- 
dicating with H the vertical scale-height of the gaseous 
component of the disk and with a the orbital distance to 
the star, hydrostatic equilibrium in the vertical direction re- 
quires that 



(3) 



in which p = p(a, z) is the gas pressure and M* is the 
mass of the star. Writing the pressure as p = c 2 p, ap- 
proximating c as being independent of height, and integrat- 
ing over the disk thickness one finds that c 2 ~ H 2 Q 2 , 
where ft = s/GM^/a 3 is the Keplerian angular veloc- 
ity of the gas. (A more accurate determination of the 
disk's rotation rate can be derived from imposing hydro- 
static equilibrium in the radial direction, which results in 
Q 2 ~ (GM*/a 3 )[l - (H/a) 2 }). Thus, equation (2) gives 



fir 



(#)ft 



(4) 



Assuming typical values at a « 5 AU for p and H/a around 
a solar-mass star of 10 -10 gcm -3 and 0.05, respectively, 
and expressing both the density and radius of the particle in 
cgs units, ClTf ~ 10 _3 p sp R sp . This estimate is applicable 
up to values of R sp less than tens of centimeters, i.e., on the 
order of the mean free path of gas molecules (oc 1/p) under 
the adopted gas conditions. 

Equation (1) has the form a damped harmonic oscillator 
and can be integrated once initial conditions are provided 
for position and velocity. Setting z = zo and z = dz / dt = 

at time zero and taking into account the inequality Qrf <C 

1 found above, the solution can be approximated as 



z~zo[e 



-tn* 



r f 



)■ 



(5) 



The second term in the above solution is a fast transient 
that rapidly decays to zero, hence we can further approx- 
imate the solution as z « zoe~ tQ T f . Therefore, in or- 
der for the altitude above the midplane to reduce by more 
than 99% of its initial value, Sit > 5/(Qt/) ~ 5x 
10 3 /(p sp i? sp ), in which again p sp and R sp are expressed 
in cgs units. Micron- sized particles would take millions of 
orbital periods (r ro t = 2tt/Q) to settle, which suggests that 
they first need to aggregate and grow into larger particles. 



Centimeter- sized icy /rocky aggregates would require only 
on the order of thousands of orbits to approach the disk's 
midplane. 

A thin layer of solid material may thus accumulate at the 
disk's midplane in a relatively short amount of time com- 
pared to disk lifetimes. If this layer is sufficiently dense, 
clumps may form through gravitational instabilities within 
the layer (e.g., Goldreich and Ward, 1973). Such process 
may produce kilometer- sized bodies, known as planetesi- 
mals. 

The assumption that the gas vertical velocity is small 
compared to the particle vertical velocity is, however, only 
valid in absence of sustained turbulent motions. In gen- 
eral, protoplanetary disks are expected to be somewhat tur- 
bulent and moderate amounts of turbulence could affect the 
settling timescales of small grains (R sp < 1cm). Ignor- 
ing the second term in equation (5), particle speeds are 
z & —zfl 2 Tf. Assuming a turbulent kinematic viscosity 
of the form v = aHc (Shakura and Syunyaev, 1973) and 
indicating with A tur b the vertical mixing length (the typical 
size of eddies), vertical gas speeds due to turbulence can be 
estimated as |i turb | ~ f/A tU rb or |i turb | ~ aHc/X turh = 
aH 2 ft I Aturb 5 and thus 



^turb 



turb 



H 2 



fir 



(6) 



Since z A tur b / H 2 < 1 (eddies cannot be larger than the 
disk's thickness), we have that |i/£turb| < Slrf/a. No- 
tice that in order for turbulent motions to be subsonic 
(i.e., |iturb| < c )> the turbulent kinematic viscosity v < 
c Aturb < cH and thus a < 1. Under the disk's condi- 
tions adopted above, the ratio of the particle's settling ve- 
locity to the gas vertical turbulent speed is then | z / i tur b | < 
10" 3 p sp Rsp/a (p sp and R sp are in cgs units). 

Therefore, if the turbulence parameter 1 a > 10 -3 , the 
settling time of centimeter- sized particles may be affected 
since vertical turbulent mixing could influence their verti- 
cal transport. Studies of grain settling in the presence of 
turbulent motions indicate that particles tend to concentrate 
in stagnant regions of the flow and that concentrations are 
size-dependent, which may lead to the accumulation of sub- 
centimeter- sized particles (Cuzzi et al., 2001). Previous 
growth of these agglomerates would still rely on sticking 
collisions of smaller particles. 

The growth from centimeter- sized to kilometer- sized 
bodies is still a poorly understood process and an ac- 
tive field of both theoretical and experimental research. 
Nonetheless, there is observational evidence that it does 
occur in nature: dust particles are observed in debris disks 
around other stars and small bodies of tens to hundreds of 
kilometers in size are observed in the Kuiper belt around 
the Sun. 



1 Based on observed gas accretion rates and other properties of young stellar 



objects with disks, typical values of a are in the range from < 
~ 0.1 (e.g., Hueso and Guillot, 2005; Isella et al, 2009). 
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2.2. From Planetesimals to Planetary Cores 

We shall now assume that the solid component of the 
protoplanetary disk has had time to agglomerate into plan- 
etesimals, rocky /icy bodies of a kilometer (or larger) in 
size. These objects are massive enough (> 10 15 -10 16 g ~ 
10 -12 M®) that they may gravitationally interact with their 
neighbors and perturb their velocities. As a result of these 
interactions, planetesimals become prone to collisions. 

Although the outcome of a collision between two plan- 
etesimals depends upon their relative velocity, we assume 
that collisions among such objects result in mergers rather 
than fragmentation. Under this assumption, the growth rate 
of a planetesimal of mass M s can be written as (Safronov, 
1969) 

dM s 2 

—jj- = ttR v re i p s F g , (7) 

where R is the planetesimal radius, v re \ is the relative ve- 
locity between the two impacting bodies, p s is the volume 
density of the solid component of the disk. The product 
v re \ p s represents a mass flux, that is the amount of solid 
material sweeping across the target planetesimal per unit 
time and unit area. The quantity F g = (R e ^/R) 2 is a grav- 
itational enhancement factor, which is the ratio of the ef- 
fective cross-section (ttR^) of the accreting planetesimal 
to its geometrical cross-section (ttR 2 ). This factor accounts 
for the ability of the growing body to bend towards itself 
the trajectories of other, sufficiently close, planetesimals. 

It is often useful to express the accretion rate dM s /dt 
in terms of the surface density of the solid material E s = 
J p s dz ~ H s p s , with H s being the vertical thickness of 
the planetesimal disk. In order to do so, one can assume 
that the gravitational force exerted by the star is the major 
component of the force acting on planetesimals in the ver- 
tical direction of the disk (as in equations (1) and (3)) and 
that relative velocities between planetesimals are isotropic. 
Thus, H s is of order v re \/Q, where Q is the Keplerian angu- 
lar velocity of the growing planetesimal. Equation (7) can 
then be cast in the form dM s /dt = ttR 2 ^ E s F g . Since 
the angular velocity along a Keplerian orbit is Q oc a -3 / 2 , 
the accretion rate of planetesimals is slower farther from 
the star (neglecting variations of T, s F g with distance). This 
would imply that dM s /dt at the current location of Uranus 
(19.2 AU) was about 7 times as small as it was at the current 
location of Jupiter (5.2 AU). 

If we neglect the collective gravitational action of the 
other planetesimals and that of the star during an encounter, 
two interacting bodies can be described in the framework 
of a two-body problem. Assuming that the target planetesi- 
mal has already grown somewhat larger than the neighbor- 
ing bodies, hence becoming a planetary embryo, the im- 
pacting body can be thought of as a projectile. We can 
therefore use the approximation that the embryo sits on 
the center of mass of the two-body system. In the rest 
frame of the embryo, conservation of the specific angular 
momentum (i.e., angular momentum per unit mass) reads 
Rett v re \ ~ R vt & , where v t & is the relative velocity for 



a tangential approach (when the projectile grazes the em- 
bryo). Conservation of specific energy during the collision 
requires that v 2 el « v 2 a — 2GM S /R and therefore the grav- 
itational enhancement factor, i.e., the ratio (R e ^/R) 2 , can 
be cast in the form 

Fg = 1 + ^ > 

V rel 

in which v esc = ^2GM s j R is the escape velocity from the 
surface of the target planetesimal. Notice that, if the plan- 
etesimal radius is not negligible compared to the embryo 
radius, the radius R in equation (7) and in v esc should be 
replaced by the sum of the two radii. 

If relative velocities are high and v re \ ^> v esc , then 
F g « 1 (i.e., i? e ff « R) and equation (7) yields dM s /dt oc 

R 2 oc Ms^ 3 . If relative velocities are low and v re \ <C v esc , 
then F g « v 2 sc /v 2 el oc R 2 and dM s /dt oc i? 4 oc M s 4/3 . If 
the escape velocity is very much larger than the relative ve- 
locity, then three-body effects (star, planetary embryo, and 
planetesimal), neglected in deriving equation (8), must be 
taken into account to compute F g (Greenzweig and Lis- 
sauer, 1992). 

In the high relative velocity regime, the growth timescale 
r s = M s (dM s / dt)~ l (i.e., the time it takes for the mass of 
the embryo to increase by a factor e « 2.7), is proportional 

1 /3 

to M s and therefore the growth rate, l/r s , of an embryo 
reduces as it grows larger. This implies an orderly growth 
of large planetesimals, which tend to attain similar masses. 
From equation (7) cast in terms of solids' surface density, 
one obtains that the timescale r og necessary to build a large 
embryo, or planetary core, of mass M can be estimated as 

To assemble a body of mass M ^ 1.6 x 10 -4 M e (about 
as massive as Ceres, the largest object in the asteroid belt) 
within a solids' disk of density 10 gem -2 , orderly growth 
would require a few 10 6 orbital periods, which already rep- 
resents a fairly long timescale compared to lifetime of pro- 
toplanetary disks around solar-mass stars (< 10 7 years). 
In the low relative velocity regime, the growth rate of 

1 /3 

an embryo is l/r s oc M s , and thus the larger the plan- 
etary embryo the faster it grows, a process known as run- 
away growth. During runaway growth, the largest embryo 
grows faster than any other embryo within its accretion re- 
gion. Although a two-body approximation (equation (8)) 
yields an unlimited gravitational enhancement factor as the 
ratio v re \/v esc — >> 0, gravitational scattering due to three- 
body effects set a limit to F g , which cannot exceed values 
much beyond several thousands (Lissauer, 1993). 

During the assembly of a planetary core through the 
growth of an embryo, relative velocities among planetesi- 
mals play a crucial role in determining the accretion rates, 
as indicated by equations (7) and (8). The velocity distribu- 
tion of a swarm of planetesimals is affected by a number of 
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physical processes, such as elastic and inelastic collisions, 
gravitational scattering, and frictional drag by the gas, the 
results of which can be highly stochastic. 

If the orbit radial excursion of planetesimals, relative to 
the orbit of an embryo, is on the order of the embryo's Hill 
radius 2 



M s 
3M* 



1/3 



(10) 



or smaller, then the accretion is said to be shear-dominated 
because growth is dictated by Keplerian shear in the disk, 
rather than by planetesimals' random velocities. If the em- 
bryo's orbit is nearly circular, this situation requires that or- 
bital eccentricities and inclinations of planetesimals should 
be ~ Rn/a, or smaller. The relative velocity between an 
embryo and a planetesimal traveling on a circular Keplerian 
orbit with radii a and a + Aa (Aa/a <C 1), respectively, is 
v re} ~ a\n(a + Aa)-Q(a)\ « aft (a) | [1 -3Aa/(2a)] - 1|, 
hence |i> re i| ~ 3ftAa/2. If we approximate Aa as the half- 
width of the region within which the gravity field of the 
embryo dominates over that of nearby embryos, the accre- 
tion rate is dM s /dt oc M^ 3 ftAa p s F g . We will see in 
Section 2.3 that Aa is generally proportional to the em- 
bryo's Hill radius, R^. Since |i> re i| ~ flRu and p s ~ 
T^s/Hs ~ E s / Ru , the planetesimal accretion rate becomes 

dM s jdt oc Ms /3 nj: s F g . 

This result is formally the same as the accretion rate in 
the orderly regime (in which v 2 sc /v 2 el <C 1 and F g w 1), 
except that now v 2 sc /v 2 el ~ 6Rn/R ^> 1, and therefore 
F g ^> 1. Accordingly, the growth timescale (in units of 
1/ft) is given by the righthand side of equation (9) di- 
vided by the gravitational enhancement factor, F g . The 
growth rate of an embryo is much larger than it is during 
the orderly growth phase (since F g ^> 1), but it reduces 
as the embryo mass grows larger (see equation (9)). This 
phase of growth, often referred to as oligarchic, may lead 
to the formation of massive embryos at regular intervals 
in semi-major axis. Notice that, neglecting variations of 
T, S F 9 with orbital distance, as in the orderly growth regime 
dM s jdt oc M^ 2 a~ 3 ^ 2 . This implies that the accretion rate 
of planetesimals reduces as the stellar mass decreases or as 
the distance from the star increases. 

2.3. Isolation Mass of Planetary Cores 

The rapid (runaway/oligarchic) growth of a planetary 
embryo continues until its neighborhood, or feeding zone, 
has been substantially cleared of planetesimals (Lissauer, 
1987). The feeding zone represents the domain within 
which the dominant embryo is able to significantly deflect 
the paths of other planetesimals towards itself. The radial 



2 The Hill radius, Rji , represents the distance of the Lagrange point L\ from 
the secondary in the circular restricted three-body problem (e.g., Murray 
and Dermott, 2000). It provides a rough measure of the distance from 
the secondary beyond which the gravitational attraction of the primary and 
centrifugal effects prevail over the gravity of the secondary. Note however 
that this region, which identifies the Roche lobe, is not a sphere and its 
volume is about a third that of a sphere of radius Rh (e.g.,Eggleton, 1983). 



extent of this region is several Hill radii, Ru, which can be 
understood by recalling the definition of Ryl . This is the dis- 
tance of the equilibrium point from mass M s (the embryo), 
on the line connecting (and in between) masses M* and M s , 
in a reference frame rotating at angular velocity ft. The 
force balance for a planetesimal of mass much smaller than 
both M* and M s requires that M S R^ 2 - M*(a - Ru)~ 2 + 
ft 2 (a — Rh) = 0. Hence, to leading order in R^/a, 
MsR^ 2 ~ M*a- 2 (1 + 2R H /a) + M*a~ 2 (l - R u /a) = 0, 
whose solution is the righthand side of equation (10). 

Once the feeding zone has been severely depleted, the 
planetary embryo becomes nearly isolated. Numerical N- 
body simulations suggest that isolation of an embryo oc- 
curs once a region of width Aa ~ bRn, where b w 4, on 
each side of the embryo's orbit becomes nearly emptied of 
planetesimals (Kokubo and Ida, 2000). The isolation mass, 
Mi SO , of a planetary embryo can therefore be calculated as 
M iso ~ 47r a Aa E s ~ 47r a b Ryl H s . Notice that here E s 
refers to the initial value of the solids' surface density. By 
using equation (10), the isolation mass can be written as 



M] f 



/ (47ra 2 6E s ) a 



3M* 



(11) 



According to equation (11), a planetary core of mass 
~ 11 M e (if 6^4) would become nearly isolated at 
5.2 AU from a solar-mass star if the local surface density 
of solids were equal to 10 gem -2 . This value for E s is a 
couple of times as large as that predicted for a minimum 
mass proto-solar nebula 3 (MMSN). Additionally, the isola- 
tion mass increases with distance from the star, although 
this increase can be somewhat compensated for by a re- 
duced surface density of solids, which is expected to de- 
crease as a increases. 

An order-of-magnitude estimate of the timescale, T{ so , 
required to reach isolation can be obtained by taking the 
ratio of equation (1 1) to the planetesimal accretion rate (ex- 
pressed in terms of S s , see Section 2.2), which yields 



ISO 

f7 



3 \ 2 / 3 



a 2 S s 1 MW 



(12) 



where C iso = (2tt) 2 / 3 (2/ v / 3) 5 / 3 ^/b/^ ~ 2.44^. Notice 
that for given values of F g and E s , this timescale increases 
with increasing distance from the star and decreases with 
increasing stellar mass: T[ so oc a 5 / 2 M^ 2 ^ 3 . At the current 
location of Jupiter around a solar-mass star, ftr{ so ~ 3 x 
10 8 /F^ , if S s is 10 g cm -2 . Assuming a situation in which 
the gravitational enhancement factor is of order 1000 over 
the entire accretion epoch, Ti SO is several tens of thousands 
of orbital periods or several 10 5 years. 



3 This is defined as the amount of heavy elements (heavier than Helium) 
observed in the planets and minor bodies of the solar system (mostly con- 
tained in giant planets) augmented by an amount of gas such to render the 
proto-solar nebula composition equal to that of the young Sun. Such def- 
inition constrains the total mass, ~ 0.02 Mq , and, to a lesser extent, the 
surface density distribution of the minimum mass solar nebula (see, e.g., 
Davis, 2005). 
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In deriving equation (12) we assumed that the surface 
density of solids, E s , is comparable to its initial value and 
the geometrical cross-section of the planet only depends on 
the core radius, R. Yet, the surface density of solids drops 
if accreted planetesimals are not replaced by others arriving 
from outside of the core's feeding zone, which would op- 
erate towards increasing T{ so . Moreover, during its growth 
toward isolation, a planetary core also accretes gas from the 
disk, although at a much smaller rate than it accretes solids. 
Once the atmospheric envelope becomes massive enough 
to dissipate the kinetic energy of incoming planetesimals 
via drag friction, the geometrical cross-section of the pro- 
toplanet in equation (7) becomes substantially larger than 
ttR 2 . This effect would operate towards reducing T[ so . 

A planetary core does not necessarily stop growing once 
attaining the mass given by equation (11). Perturbations 
among planetesimals and other embryos can supply addi- 
tional solid material to the core's feeding zone via scatter- 
ing. Furthermore, a planet in excess of a Mars mass can 
exert gravitational torques on the surrounding gas. The disk 
responds by exerting the same amount of torques onto the 
planet, which modify the planet's orbital angular momen- 
tum, forcing it to radially migrate within the disk. As a 
result of this radial displacement, a planet may reach disk 
regions that still contain planetesimals. However, as the 
planet's mass grows, so does its ability to scatter planetesi- 
mals away from the orbital path or to trap them into mean- 
motion resonances (so that the ratio of a planetesimal's or- 
bital period to the planet's orbital period is a rational num- 
ber). 

The order-of-magnitude estimates given above for the 
isolation mass and the timescale to reach it neglect many 
aspects of the physical processes involved in the growth of 
a planetary core. Nonetheless, detailed calculations of gi- 
ant planet formation by core accretion and gas capture indi- 
cate that those estimates are valid under appropriate con- 
ditions. Calculations that start from a planetary embryo 
of M s ~ 0.1 M e (about equal to the mass of Mars) or- 
biting a solar-mass star at 5.2 AU and undergoing rapid 
growth within a planetesimal disk of initial surface den- 
sity E s « 10 gem -2 , show that a planetary core becomes 
nearly isolated within less than half a million years when 
Mi SO ^11 M e (Pollack et al., 1996), as also indicated by 
the solid line in left panel of Figure 1 (see Figure's caption 
for further details). 

2 .4 . Growth of Thermally-Regulated Envelopes 

Gas can accrete onto a planetary embryo when the ther- 
mal energy is smaller than the gravitational energy bind- 
ing the gas to the embryo. This condition is satisfied 
when the escape velocity from the surface of the embryo, 
v esc (see Section 2.2), exceeds the local thermal speed 
of the disk's gas (^/S/ttc), which occurs when M s > 
4(H/a) 2 M*R/(7ra),or 



At distances of several AU from a solar-mass star, relatively 
small bodies (M s ~ 0.01 M®) can retain an atmosphere. In 
these early phases, the atmospheric gas is optically thin and 
thermal energy released by impacting planetesimals can be 
readily radiated away, allowing for contraction of the enve- 
lope and for additional gas to be accreted. 

Prior to achieving isolation from the planetesimal disk, 
the accretion rate of solids is much larger than that of gas. 
The growing core mass forces the envelope to contract and 
more gas can be collected from the surrounding gaseous 
disk. However, as the envelope grows more massive, it 
becomes increasingly optically thick to its own radiation, 
which therefore cannot escape to outer space as easily as in 
earlier phases. As a consequence, envelope's temperature 
and density exceed those of the local disk's gas. 

The pressure gradient that builds up in the envelope 
opposes gravitational contraction, preventing accretion of 
large amounts of gas. On the one hand, ongoing accretion 
of solid material and growth of the core help envelope con- 
traction, on the other the gravitational energy released in the 
envelope by gas compression and that supplied by accreted 
planetesimals help maintain a relatively large pressure gra- 
dient. Therefore, contraction becomes self-regulated. Once 
the planet achieves (near) isolation from planetesimals, the 
accretion rate of solids becomes small and accretion of gas 
continues to the extent allowed by envelope compression. 

The envelope enters a stage of quasistatic contraction 
that can be characterized by long evolution timescales. If 
the accretion of solids and gas was negligibly small, the 
timescale over which the envelope contracts would be re- 
lated to the ratio of the gravitational energy released by con- 
traction, |£^g r av | , to the envelope's luminosity, L. This ratio 
defines the Kelvin-Helmholtz timescale tkh = \Eg rdLV \/L. 
Indicating the protoplanet's core and envelope masses with 
M c and M e , respectively, |£ grav | ~ GM c M e /R, R being 
the average radius comprising most the protoplanet's mass, 
and thus the contraction timescale is 



tkh 



GM c M e 
RL 



(14) 



M, > 



M? ( H 



(13) 



Using values from Figures 1 and 2 around the middle of the 
planetesimal isolation phase, tkh would be of order 10 5 
years. However, for the case illustrated in the Figures, most 
of the luminosity produced during this slow growth phase 
is due to gravitational energy generated by accretion rather 
than by contraction of the envelope. 

The length of this epoch depends on several factors, but 
principally on the opacity of the envelope and on the solids' 
accretion rate. The more opaque the gas is to outgoing ra- 
diation, the less able the envelope is to dissipate the energy 
provided by gas compression and/or by continued planetes- 
imal accretion. Moreover, accreted solids can be consumed 
by ablation during their atmospheric entry, releasing dust 
and increasing the envelope's opacity. Therefore, a reduced 
accretion of solids shortens this epoch. As shown in the left 
panel of Figure 1 , the slow contraction phase, from the time 
the core mass reaches near-isolation from surrounding plan- 
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Figure 1. Formation models of a giant planet by core nucleated accretion need to take into account many physical effects. These 
include: i) Calculation of the planetesimal accretion rate onto the protoplanet; u) Calculation of the interaction via gas drag of impacting 
planetesimals with the protoplanet's gaseous envelope; in) Thermodynamics calculation of the envelope structure; iv) Calculation of 
the gas accretion rate from the protoplanetary disk onto the protoplanet during the phase of slow envelope growth; v) Hydrodynamics 
calculation of constraints on envelope size and disk-limited gas accretion rates during the phase of runaway gas accretion. Some results 
from one such calculation are reported here for a model of Jupiter's formation (Lissauer et al., 2009). This model assumes that the 
planetesimal disk has S s « 10 g cm -2 and the grain opacity in the protoplanet's envelope is 2% the value of the interstellar medium. 
The gaseous disk is assumed to dissipate within 3 Myr. The left panel shows the mass of solids (solid line), of gas (dotted line), and the 
total mass of the planet (dot-dashed line) as functions of time. The right panel displays the radius of the planet (solid line) and that of 
the planet's solid core (dashed line) from the same model. 



10" 3 F 




12 3 

t (Myr) 

Figure 2. Luminosity as a function of time from the model re- 
ported in Figure 1. From the time of the first luminosity peak 
(t ~ 0.35 Myr) to the middle of the planetesimal isolation phase 
(t ~ 1.6 Myr), the luminosity of the protoplanet decreases by 
more than a factor of 25. This is mainly due to the reduced accre- 
tion rate of solids. 

etesimals to the time the rapid gas accretion phase begins, 
lasts about 2 Myr. The assumption made for this particular 
model is that the grain opacity in the protoplanet's enve- 
lope is 2% that of the interstellar medium and that plan- 
etesimals continue to accrete throughout the slow contrac- 
tion phase. All else being equal, this time-span increases to 
6 Myr when the full interstellar opacity is adopted and de- 
creases to « 1 Myr if planetesimal accretion ceases entirely 



once the isolation mass is achieved (Hubickyj et al., 2005). 

It thus appears that the duration of the slow contrac- 
tion phase represents the main uncertainty in determining 
whether a giant planet can form before the gaseous com- 
ponent of the disk dissipates (i.e., within a few to several 
millions of years; see the chapter by Roberge and Kamp). 
However, it offers a natural explanation for the formation 
of planets such as Uranus and Neptune in our solar sys- 
tem, which do not posses very massive envelopes. These 
planets may have not been able, because of core and enve- 
lope conditions, or may have not had enough time, because 
of a dissipating disk, to acquire a large envelope. In addi- 
tion, static atmospheric models indicate that the removal of 
grains by growth and settling from the radiative zone of a 
planet's envelope may significantly reduce their contribu- 
tion to the opacity (Movshovitz and Podolak, 2008), hence 
shortening the contraction timescale of the envelope. 

Quasi- static calculations of the envelope thermal struc- 
ture suggest that the growth timescale of a protoplanet's en- 
velope r e = M(dM e /dt)~ 1 can be cast in the following 
parametric form (e.g., Ida and Lin, 2008) 

in which M = M c + M e is the total planet's mass and 
M c > Mi SQ . The timescale f e depends on several factors, 
but mostly on the opacity of the envelope. For the model 
reported in Figure 1, f e ~ 10 10 years and £ ~ 3. No- 
tice, though, that many factors affect the values of param- 
eters £ and, especially, of f e . In general, applicability of 
equation (15) requires some previous knowledge of the en- 
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velope thermal conditions, for example through theoretical 
arguments or numerical modeling. When M e > M c , the 
inner (densest) parts of the envelope, which contain most of 
the mass, effectively contribute to the compression of the 
outer parts. Yet, around this stage, the envelope may be- 
come gravitationally unstable and collapse (see Figure 1). 

2.5. Critical Core Mass and Envelope Collapse 

In order to gain some insight into the conditions lead- 
ing to envelope collapse, we can construct simplified ana- 
lytical models that qualitatively describe the physical state 
of the envelope. The first assumption we shall make is 
that a slowly contracting envelope evolves through stages 
of quasihydrostatic equilibria. Early studies of a static enve- 
lope surrounding a planetary core (e.g.,Perri and Cameron, 
191 A) suggested that there may not be any static structure 
beyond some value of the core mass. Past such a critical 
core mass, the extended outer portions of the gaseous enve- 
lope become hydrodynamically unstable and collapse onto 
the core. Simplified models of hydrostatic envelopes can 
shed some light over the concept of critical core mass. One 
such model was used by Stevenson (1982) to interpret re- 
sults of static envelope calculations performed by Mizuno 
(1980). 

Let us assume that a gaseous envelope of mass M e sur- 
rounds a core of mass M c and that the envelope is spheri- 
cally symmetric around the core. Hydrostatic equilibrium 
in an envelope shell of radius r and thickness dr implies 
that there is a balance between gravitational force and pres- 
sure acting at the shell surface boundaries. Therefore, the 
condition for hydrostatic equilibrium reads 



dp 
dr 



Gmp 



(16) 



where p and p are, respectively, the pressure and density at 
radius r, and m is the mass within r. 

An important assumption we shall make is that energy 
is transported through the envelope via radiative diffusion 
(i.e., this model applies to a radiative envelope) and that 
the envelope is stable against convection. In general, this 
may be true only in the outer layers of the envelope (Bo- 
denheimer and Pollack, 1986). 

The amount of radiation energy transported through an 
envelope shell per unit surface area and unit time, i.e., the 
radiation flux, is F ra d = —VdE^jdr, where the radia- 
tion energy density i? ra d oc T 4 and T is the temperature. 
The radiation diffusion coefficient, V, is proportional to the 
speed of light times the mean free path of photons, 1 / (hp) , 
where R is a frequency-integrated opacity (or mass absorp- 
tion coefficient). Hence, for an envelope shell of mass 
dm, the radiation flux is also equal to (e.g., Kippenhahn 
andWeigert, 1990) F rad = -167rr 2 cr rad / (3R)(dT 4 /dm), 
where a ra d is the Stefan-Boltzmann constant. If the lumi- 
nosity through the envelope, L, is nearly constant then one 
can write F ra d = 47rr 2 L, and thus 

dT 4 3RL 
dm 647T 2 cr ra d r 4 ' 



(17) 



Another assumption we shall rely upon is that the gra- 
dient V = d In T/d In p is nearly constant in the envelope. 
The temperature gradient can thus be cast in the form 



dT 
dr 



V. 



(18) 



In order for energy to be transported via radiation alone, 
the gradient V must be smaller than its adiabatic value, 
V a d • For a perfect monoatomic gas and negligible radia- 
tion pressure, V a d = 2/5. If V > V a d, the envelope be- 
comes convectively unstable. By substituting equation (16) 
in place of the pressure gradient in equation (18) and using 
the equation of state for a perfect gas p/p = 1ZT / p (JZ is 
the gas constant), we have that dT/dr ~ -pGMV / (Jlr 2 ), 
where M = M c + M e , and hence T - pGMV/(Kr). In 
the above relations, p indicates the mean-molecular weight, 
that is the average number of atomic mass units per gas par- 
ticle. Notice that the approximation m ~ M made above 
determines how deep in the envelope this temperature strat- 
ification is applicable. 

Equation (17) can be written as 



dm 
dr 



4 4 7T^ CT rad 

3 RL 



4 rri [ 

r T 



dT 
dr 



from which it follows that 



dm _ 4 4 tt 2 a rad / pGMV Y 
~dr 3RL I ^ J 



(19) 



(20) 



Integrating both sides of this equation along r , from the core 
radius, R c , to the envelope radius, R, and neglecting opac- 
ity variations 4 , one obtains 



M-M C = A ^-^M\ 
n L 



(21) 



in which A = (n 2 a rad / 3) (4G/1Z) 4 In (R / R c ) . 

We seek for the largest value of the core mass that can 
have a fully radiative envelope in hydrostatic equilibrium. 
This can be obtained from equation (21), which gives an 
explicit function M c = M C (M) , if L is independent of both 
M and M c . The value of M for which M c is maximum is 
given by 

1/3 



1 RL 
4^ (/iV) 4 



(22) 



Physically significant solutions of equation (21) can only 
be obtained for M/M cr < 4 1 / 3 , and therefore the mass 
of strictly hydrostatic envelopes does not exceed 4 1 / 3 M cr . 
The critical core mass for radiative envelopes is found by 
substituting equation (22) in equation (21), which gives 



(23) 



4 A more general opacity law, given by the product of some power of p and 
some power of T, would result in the same relation as equation (21), but 
with a modified form of coefficient Ao (Wuchterl et al., 2000). 
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According to these simple analytical arguments, the crit- 
ical core mass, M£ r , is independent of the local density and 
temperature of the protoplanetary disk, it is weakly depen- 
dent on envelope opacity and luminosity, but more strongly 
dependent on /i and V. Using radii and luminosity values 
before collapse from Figures 1 and 2, \i = 2.25, V = 0.235 
(Kippenhahn and Weigert, 1990), and R = 0.02 cm 2 g -1 , 
one obtains M£ r ~ 16 M e and a total mass of ~ 21 M 
(recall that these values would apply to a strictly hydrostatic 
and fully radiative envelope). 

As mentioned in Section 2.4, most of the protoplanet's 
luminosity during the slow growth phase may be due to 
gravitational energy released by accreted planetesimals. In 
these circumstances the luminosity can be approximated as 
L ~ (GM c /R c )(dM c /dt), in which dM c /dt is the plan- 
etesimal accretion rate, thus L oc M^ 3 (dM c /dt). Sub- 
stituting in equation (21), one obtains M c as an implicit 
function of M. By differentiation of this function, one can 
readily show that again M c cr = (3/4) M cr . 

The deepest and densest layers of an envelope are typ- 
ically convective (e.g., Bodenheimer and Pollack, 1986). 
Simple analytical models of hydrostatic envelopes can also 
be constructed under the hypothesis that energy is trans- 
ported only via convection, as that proposed by Wuchterl 
(1993). By employing the adiabatic gradient V a d, the tem- 
perature gradient becomes dT/dr ~ -/iGMV a d/(^ 2 ). 
Instead of an energy equation, in a fully convective enve- 
lope one can use the polytropic law T/Td = (p/pd) f , 
in which Tjj and pD indicate values where the envelope 
merges with the disk. In a perfect monoatomic gas with 
negligible radiation pressure, T = (5/3) V a d- By differen- 
tiation of the polytropic law, it follows that 



T = r (-)(- 

dp \pDj \PD 



(r-i) 



(24) 



which combined with the temperature gradient above yields 



Pd ) \Pd ) dr 



f pGMV ad 
V 1Zr 2 



(25) 



Equation (25) can be integrated in radius to obtain the 
envelope's density stratification 



The approximation m ~ M, introduced above to derive the 
temperature gradient, again determines the envelope depths 
to which this density stratification applies. In equation (26) 
the constant of integration is set to zero by choosing as outer 
envelope radius R = /iGMV a d/(KT])). By integrating 
the quantity Airr 2 pdr from the core radius, R c , to R we get 



M-M c = B p D 



3T-1 



^Vad 



M 6 



(27) 



where B = ^{G/llf [l - (R c /R) {3 ~ 1/r) ] , which can 
be approximated to 4tt(G/1Z) 3 since R c <C R and 3 — 
1/T > 1 when radiation pressure can be ignored. 



The largest core mass that bears a fully convective enve- 
lope in hydrostatic equilibrium corresponds to a total mass 



1 /3T- 1\ / n 



12tt 



U^Vad 



_ D 

PD 



1/2 



(28) 



In order for equation (27) to admit physically significant 
solutions, M/M cr < \/3. According to these simple ar- 
guments, strictly hydrostatic and convective envelopes are 
limited in mass to \/3M cr . From equations (27) and (28), 
one finds that the critical core mass of a convective envelope 
is 

M c cr = -M cr . (29) 
o 

Unlike the solution for radiative envelopes, the critical 
core mass for convective envelopes explicitly depends on 
the disk density and temperature. Using values pD ~ 
10 -10 gcm -3 and Tjj ~ 100 K, which would apply to a 
protoplanetary disk at ~ 5 AU from a solar-mass star if 
there was no planet, equation (28) gives M cr ~ 170 M e . 
However, once M achieves several tens of Earth masses, 
both po and Tp are affected by disk-planet interactions 
and hence depend on M. In these cases, because of the 
gravitational perturbation produced by the protoplanet, lo- 
cal disk densities as well as temperatures can be smaller 
than the corresponding values in absence of the protoplanet 
(D'Angelo et al., 2003a), if effects of stellar irradiation can 
be neglected. 

The two analytic solutions for fully radiative and fully 
convective envelopes can be combined to obtain an estimate 
of the critical core mass for composite convective-radiative 
envelopes. Following the same line of argument, we con- 
sider a two-layer model in which the "convective" solution 
applies to the inner convective layer that extends from the 
core radius, R c , to the radius i? c ib> which plays the role of 
envelope radius (R) in the convective solution given above. 
Accordingly, temperature T c ib and density p c ib at the con- 
vective layer boundary, i? c ib> replace the disk values T& 
and po, respectively. Then, indicating with M c ib the mass 
within Rdb (core mass plus mass of the convective layer), 
equation (29) gives a critical value M™ h = (3/2)M£ r . 

The "radiative" solution applies to the outer radiative 
layer, whose boundaries are at the convective layer (outer) 
radius i? c it>, which plays the role of core radius (R c ) in 
the radiative solution above, and at the envelope radius 
R. Thus, the mass within the convective layer boundary, 
M c ib, replaces the core mass (M c ) in the radiative solu- 
tion above. Equation (23) gives the critical value of the 
total mass M cr = (4/3)M£ r b . Therefore, the largest planet 
mass for which both the convective and radiative layers, and 
hence the entire envelope, can be in hydrostatic equilibrium 
is 

M cr = \ M c Y b = 2M c cr . (30) 

According to equation (30), a composite convective- 
radiative envelope collapses once the envelope mass is 
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equal to the core mass. This prediction of the simple an- 
alytic model is in agreement with the results from the de- 
tailed calculation shown in the left panel of Figure 1 . Equa- 
tion (22) can still be used to estimate M cr , where now 
Ao oc In (R/R c ih). At the epoch of envelope collapse, 
the calculation in Figure 1 gives a ratio R/Rd^ w 4.74 
which, together with the other values adopted above, yields 
M cr ~ 32 M e and M c cr ~ 16 M e , consistent with the re- 
sults illustrated in the Figure's left panel. 

2.6. Disk-Limited Gas Accretion Rates 

During the collapse phase, large amounts of gas can be 
accreted from the protoplanetary disk. The rate at which 
the protoplanet's envelope grows soon becomes very large. 
When the total mass is M ~ 20 M®, equation (15) yields 
an envelope growth timescale of order 10 6 years. But 
once the envelope mass, M e , exceeds the core mass (equa- 
tion (30)), this length of time can shorten by orders of 
magnitude. In the model illustrated in Figure 1, when 
M e ~ 2 M c the growth timescale becomes of order 10 4 
years. 

The protoplanetary disk feeds gas to the planet's vicin- 
ity. Yet, the rate at which such gas supply occurs is limited. 
In this stage, since a protoplanet can basically accrete gas 
at arbitrarily large rates, gas accretion is governed by hy- 
drodynamical factors that involve tidal (i.e., gravitational) 
interactions between the growing planet and the disk . To 
derive simple estimates of disk-limited gas accretion rates 5 , 
dM/dt, we can assume that a gas parcel orbits about the 
star, at a distance a + Aa, with an orbital velocity relative 
to the planet of order QAa (see Section 2.2), where Q is 
the angular velocity of the planet and a its orbital radius. 
The protoplanet is able to capture gas within an effective 
radius R gc (i.e., 7rR^ c is its effective cross-section for gas 
capture) smaller than the disk's local thickness, H. If p 
is the gas volume density, then the mass flux through the 
planet's cross-section is pQR gc . The gas accretion rate is 
given by the mass flux times the planet's cross-section, that 

in which the volume density is expressed in terms of surface 
density (p ~ D/iT). 

There are two relevant length- scales over which a pro- 
toplanet can attract gas: the Hill and Bondi radii. The Hill 
radius, Rb, is defined in equation (10) (in which M s is re- 
placed with the total mass M) and represents a measure of 
the distance past which centrifugal forces and gravitational 
forces by the star dominate over the gravitational force ex- 
erted on the gas by the planet. The Bondi radius is 



Rb 



GM 



(32) 



5 In this phase, we designate the accretion rate of gas simply as dM/dt since 
total mass variations are overwhelmingly determined by accretion of gas, 
regardless of whether there is continued accretion of solids. 



(c = HQ is the gas sound speed) and marks the distance 
beyond which the thermal energy of the gas is larger than 
the gravitational energy binding that gas to the planet. In 
the case that Rb < Rb, forces due to pressure gradients 
within the Hill sphere prevent gas from becoming bound 
to the protoplanet. Then, the Bondi radius is the relevant 
distance over which gas can be accreted. We will refer to 
this regime as Bondi-type gas accretion. The Hill radius be- 
comes the relevant distance for gas accretion if Rb < Rb , 
which will be designated as Hill-type gas accretion regime. 
Therefore, in the general case, the effective radius for gas 
capture of a protoplanet embedded in a gaseous disk is 



R gc = min(i? B ,#H) • 



(33) 



By substituting equation (32) in equation (31), one finds 
that the protoplanet's growth rate, (dM/dt)/M, in the 
Bondi-type gas accretion regime is 

whereas, in the Hill-type gas accretion regime, the growth 
rate is found by using equation (10) (with M in place of 
M s ) and equation (31), resulting in 



1 

m 



l 



CbQ 



q 2 £ 



(#)■ 



(35) 



The dimensionless constants Cb and Cb account for small 
corrections due to effects neglected in deriving this simple 
model and should be of order unity. They can be obtained 
from direct numerical simulations. Notice that dM/dt oc 
M 3 for Bondi-type gas accretion (implying a faster than 
exponential growth) while dM/dt oc M for Hill-type gas 
accretion (implying an exponential growth). 

Bondi-type gas accretion is a very steep function of the 
disk thickness, i.e., of the local gas sound speed. If the disk 
is locally very warm, high thermal energy prevents gas from 
becoming bound and hence accreted (the planet's cross- 
section for gas accretion decreases for increasing temper- 
ature). Hill-type gas accretion is less sensitive on disk tem- 
perature because in this regime gravitational factors domi- 
nate over thermal ones. 

In general, the mass growth rate of the protoplanet, 
1/td , which is limited by disk hydrodynamics, is given by 



- = ( 1 



/TB 
/TH 



for M < M tr 
for M > M tr 



(36) 



where the transition between the two regimes occurs at the 
transition mass M tr , for which tr = tb , and thus 




(37) 



This simple gas accretion model can be compared 
against results from three-dimensional hydrodynamical 
simulations of a planet gravitationally interacting with a 
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Figure 3. Left: Mass growth rate, 1/td — (dM/dt)/M, versus mass of a protoplanet accreting gas at a limiting rate provided 
by a protoplanetary disk. These results were obtained from a three-dimensional, high-resolution hydrodynamical model of a growing 
protoplanet interacting with a disk of initial surface density at a — 5.2 AU of 100 g cm -2 , relative thickness H/ a — 0.05, and turbulence 
parameter a = 4xl0 -3 (D'Angelo and Lubow, 2008). The timescale m is in units of orbital periods. The dashed line represents 
the growth rate given by equation (36): the slanted portion corresponds to the rate in the Bondi-type gas accretion regime (1/tb with 
Cb — 2.6) and the horizontal portion corresponds to the Hill-type gas accretion rate (1/th with Cu ~ 0.9). At large planet masses, 
the growth rates drop due to the formation of the tidally produced gap in the local density distribution. Right: Average surface density 
distribution of the disk's gas, E, near the planet relative to the local unperturbed (i.e., without a planet) value, Xo, as a function of the 
planet's mass. 



disk and growing in mass at a disk-limited gas accretion 
rate. Figure 3 (left panel) shows numerical results (solid 
line) and the growth rate given by equation (36) (dashed 
line). Agreement is found for values of the constants 
Cb — 2.6 and Ch — 0.9. Note that the mass scaling, 
i.e., the slopes of the two dashed line segments, is correctly 
predicted by the simple model and so is the transition mass 
between the two regimes (M tr w 4 x 10" 5 M* , or & 14 M e 
if M* = 1 M ). The disk thickness in the example shown 
in Figure 3 is H/a = 0.05, which corresponds, for typical 
disk properties, to a local temperature T ~ 100 K. Since 
the gas sound speed c oc VT, a factor of 2 increase in 
the local temperature would reduce the planet's growth rate 
by an order of magnitude during the Bondi-type accretion 
regime. 

As the planet's mass increases, the effective radius for 
gas capture may exceed the disk thickness. Moreover, den- 
sity perturbations due to tidal interactions of the protoplanet 
with the disk are no longer negligible. The gas density 
along the planet's orbit starts to be affected (see right panel 
of Figure 3), and the simple accretion model becomes in- 
applicable beyond some value of the protoplanet mass. In 
the example shown in Figure 3, this occurs for masses 
M > 10 -4 M* , when disk-planet interactions have changed 
the local surface density by more than 20% of its unper- 
turbed value (i.e., when there is no planet in the disk), as 
illustrated in the right panel. 

2.7. Gap Formation 

A protoplanet embedded in a disk exerts a gravitational 
torque on the gas interior of a— | Aa | and exterior of a+ 1 Aa | 



whose magnitude is of order 




which leads to exchange of orbital angular momentum be- 
tween the planet and the disk (see the chapter by Lubow 
and Ida for a derivation of this torque). The factor / is typ- 
ically of order unity. It can be shown that | Aa\ is the larger 
between the Hill radius, Rn, and the disk scale-height, H. 
The sign of T 9 is positive for material lying outside of the 
planet's orbit (Aa > 0) and negative for material inside 
of it (Aa < 0). Hence, material orbiting outside of the 
planet's orbit gains angular momentum, moving towards 
larger radii, whereas material inside of the planet's orbit 
loses angular momentum, moving towards smaller radii. 
This process tends to deplete the disk of gas along the 
planet's orbit, thus forming an annular gap in the local den- 
sity distribution (e.g., Lin and Papaloizou, 1986). However, 
in a viscous (Keplerian) disk, because of differential rota- 
tion gas is also subject to a viscous torque generated by 
viscous friction between adjacent disk rings. The viscous 
torque exerted by material inside of the orbital radius a on 
material outside of a is (Lynden-B ell and Pr ingle , 1974) 

T v ^37nsa 2 Ztt. (39) 

This torque tends to smooth out density gradients and redis- 
tribute material across the planet's orbit, filling in the gap. 

A condition for gap formation requires a net gain of or- 
bital angular momentum for material outside the radius a, 
that is T 9 — T v > 0, and a net loss of orbital angular mo- 
mentum for material inside a,T g + % < 0. If H > Ru, 
the condition \T g \ > % translates into a condition for the 
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minimum planet's mass necessary to open up a density gap 
in a disk 

where the turbulence parameter a is defined in Section 2.1 . 
In case of very cold disks or very massive planets Rn > H 
and the torque inequality gives (M/M*) > irfa (H/a) 2 . 
For conditions adopted in the simulations reported in Fig- 
ure 3, substantial gas depletion is expected for M > 2 x 
10 -4 M* (or M > 60 M e for a disk around a solar-mass 
star). In accord with this estimate, the right panel of the 
Figure indicates that the average surface density along the 
planet's orbit is reduced by 40% of its unperturbed value 
when M > 2 x 10~ 4 M*. 

The timescale for gap formation can be estimated from 
the equation of motion of a thin viscous Keplerian disk, 
which also evolves under the action of gravitational torques 
exerted by an embedded object (Lin and Papaloizou, 1986) 



M/Mj 

0.01 0.10 1.00 10.00 




_ _d_ 

dt da 



aQ da v 



(41) 



M/M„ 

Figure 4. Disk-limited gas accretion rates as a function of the 
planet mass obtained from three-dimensional, high-resolution hy- 
drodynamical calculations of a planet interacting with a proto- 
planetary disk. Accretion rates are in units of the unperturbed 
surface density at the planet's orbital radius, a, and the planet's 
orbital period, r ro t . The top axis uses M* = 1 M . Filled cir- 
cles correspond to results for a disk with a turbulence parameter 
a = 4xl0 -3 (see Section 2.1). Empty circles are for a disk with 
a — 4x 10 -4 . The disk's aspect ratio is H/a — 0.05. The curves 
represent fits to the data. 



orbiting a solar-mass star at ~ 5AU in a disk, whose 
gas surface density just outside the gap region is of order 



it follows that 



gap 



V M 



lAal 



(42) 



where radius a is now interpreted as a variable indicat- 
ing the distance from the star. Equation (41) can be de- 
rived from imposing conservation of mass and angular mo- 
mentum in an axisymmetric and flat disk. Ignoring vis- 
cous torques (% < \T g \) and taking the approximation 
8/ da — >> 1/ Aa in equation (41), defining the gap formation 
timescale as r gap = ^d^/dt^ 1 and using equation (38), 100 gem" ^, may accrete gas at rates - 10- 5 -10- 2 M e per 

year, thereby reaching a Jupiter's mass in 10 4 -10 5 years. 

2.8. Final Masses of Giant Planets 

It is not yet entirely clear what all the factors are that 
determine the final mass of a giant planet. Yet, tidal trunca- 
tion of the disk by gravitational interactions with the planet 
is likely one of the main factors. As indicated in the right 
panel of Figure 3 , gaps can become quite deep around a 
Jupiter-mass planet, even at moderate values of the turbu- 
lence parameter a. However, gas can filter through the tidal 
barrier (i.e., the gap) and continue to accrete towards the 
planet. 

This can be seen in Figure 5 , which shows the gap in the 
surface density distribution around a giant planet together 
with trajectories of material moving towards the planet 
along the inner and outer edges of the gap. The trajecto- 
ries are drawn in a frame co-moving with the planet. This 
material becomes trapped in the inner parts of the planet's 
Roche lobe and is eventually accreted. As mentioned in the 
previous section, accretion through a gap can be quite ef- 
ficient. Assuming a gas surface density of ~ 100 gem -2 
just outside the gap region at ~ 5 AU, the growth timescale 
(due to disk-limited gas accretion), td = M(dM/dt)~ 1 , 
of a Jupiter-mass planet in a disk with a ~ 10 -3 and 
H/a « 0.05 would be several thousands of orbital peri- 
ods (see Figure 4). Around a solar-mass star, the tidal bar- 
rier alone would constrain the planet's final mass to about 
6-7 Mj. 



with | Aa\ equal to the larger of H and Ru- This timescale 
is typically short, on the order of tens of orbital periods. 

Once gap formation starts, accretion rates drop because 
of gas depletion around the planet. According to con- 
dition (40), lower disk viscosity allows for gas depletion 
around smaller mass planets. By using equation (37), one 
finds that for 



< 



— 



(43) 



9tt/C b \a 

gas depletion begins before the transition to a Hill-type gas 
accretion and the growth rate starts to decline for smaller 
mass planets. In a disk of thickness H/a « 0.05, this may 
occur for a less than few times 10 -4 . 

Effects of gap formation can be seen in Figure 4. The 
Figure shows disk-limited gas accretion rates, derived from 
hydrodynamical simulations (see D'Angelo et al., 2003b, 
for details), for protoplanets embedded in a moderately 
viscous (a = 4 x 10 -3 ) and in a low-viscosity (a = 
4 x 10 -4 ) disk. Notice, though, that if r gap > th, i.e., 
/C H (a 2 £/M*) > (M/M*) 2 (a/H)\ Hill-type gas accre- 
tion may persist until Rn ~ H . 

Disk-limited gas accretion rates can be quite large, un- 
less the disk's gas density is low. A Saturn-mass planet 



D'Angelo, Durisen, & Lissauer 



12 of 28 



Giant Planet Formation 



2.9 Orbital Migration of Protoplanets 



2 CORE ACCRETION MODELS 



2.1E + 00 1.4E + 02 2.8E + 02 4.2E + 02 5.6E + 02 




-2-10 1 2 

X/o. 



Figure 5. Surface density of a gaseous disk containing a pro- 
toplanet whose mass is M — 10 -3 M* (i.e., one Jupiter's mass 
for a solar-mass star). The relative disk thickness is H/a = 0.05 
and the turbulence parameter is a — 4 x 10 -3 . The planet is lo- 
cated at ( — 1, 0), the star at (0, 0), and the disk is rotating in the 
counterclockwise direction. The grey scale bar is in cgs units. The 
plot shows the density gap along the planet's orbit and the wave 
pattern generated by disk-planet interactions. The two black lines 
represent trajectories of gas in the co-moving frame of the planet. 
These gas parcels move along the inner and outer gap edges, be- 
come gravitationally bound to the planet, and are eventually ac- 
creted. 

A lower disk temperature, i.e., a smaller disk thickness, 
can increase the planet's tidal barrier and reduce gas ac- 
cretion. However, it is unclear whether accretion can be 
stopped around M ~ 1 Mj for reasonable values of H/a 
(Lubow and DAngelo, 2006). Lower disk viscosity also 
helps reduce disk-limited gas accretion rates, as indicated 
in Figure 4. Assuming a disk around a solar-mass star with 
a ~ 10 -4 , for the same example made above, the tidal bar- 
rier would limit the planet's final mass to ~ 2 Mj. 

The considerations above do not take into account the 
fact the a protoplanetary disk evaporates due to irradia- 
tion from the central star and, possibly, from other exter- 
nal sources (Hollenbach et al., 2000). Both observations 
and theory suggest that photoevaporation timescales are on 
the order of a few million years. Therefore, disk consump- 
tion by photoevaporation may become a relevant process in 
the final stages of a giant planet's growth. It is likely that 
a combination of these three factors, i.e., a reducing disk 
thickness H/a, a decaying turbulence a, and gas depletion 
by photoevaporation, all of which may take place as a proto- 
planetary disk ages, plays a significant role in determining 
final masses of gas giant planets. 

Once gas accretion stops, a giant planet evolves at a con- 
stant mass. The planet gradually cools down and contracts 



on a Kelvin-Helmholtz timescale (equation (14)). The lu- 
minosity is mostly due to gravitational energy released in 
the contraction. 

2.9. Orbital Migration of Protoplanets 

Gravitational interactions between a protoplanetary disk 
and an embedded planet result in an exchange of orbital an- 
gular momentum (see the chapter by Lubow and Ida for a 
thorough discussion on this topic). Assuming that H > Ru 
and that condition (40) for gap formation is not satisfied, 
the disk's surface density is largely unperturbed and equa- 
tion (38) applies. 

If the planet exerts a torque T g on the gas outside its or- 
bit, then conservation of angular momentum dictates that 
this gas exerts the same torque (with an opposite sign) on 
the planet. Therefore, the planet loses orbital angular mo- 
mentum to material exterior of its orbit. By the same prin- 
ciple, the planet gains orbital angular momentum from ma- 
terial interior of its orbit. The torque expression in equa- 
tion (38) is symmetrical with respect to the planet's radial 
position, which would result in a zero net torque acting on 
the planet. However, one can show that, because of global 
variations of the disk properties across the planet's orbital 
radius a, the net torque exerted on the planet by disk mate- 
rial is of order \T P \ ~ (H/a) \T g \ (Goldreich and Tremaine, 
1980; Ward, 1997), hence 

Quantity C\ depends on the radial gradients of surface 
density, temperature, and pressure of the disk across the 
planet's orbit. In typical circumstances, this quantity is neg- 
ative and of order unity (Tanaka et al., 2002). 

The orbital angular momentum of a planet on a circu- 
lar orbit is L p = Ma 2 ft. If the planet is acted upon by a 
torque T p , conservation of angular momentum imposes that 
dL p /dt — T p = 0, and thus 

1 dM lda = T p 
M dt 2a dt Ma 2 ft ' K } 

If the first term on the lefthand side of equation (45) is neg- 
ligible compared to the second, then da/dt = 2T P / (MaQ) . 
This phenomenon is referred to as orbital migration. Since 
the torque experienced by the planet is generally negative, 
the orbit shrinks. When the torque acting on a (non-gap 
opening) planet is of the type in equation (44), migration 
is said of type I. The migration timescale, a\da/dt\~ l , in 
these circumstances is 

(£)(!)'■ «« 

Orbital decay due to type I migration can be quite rapid. 
A 10 M e (M = 3 x 10" 5 M*) planetary core interacting 
with a disk, whose surface density is E ~ 100 gem -2 and 
H/a ~ 0.05, would migrate on a characteristic timescale 
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of ~ 4 x 10 4 orbital periods. For a planet orbiting at ~ 
5 AU around a solar-mass star this period would amount 
to several times 10 5 years, a factor of a few shorter than the 
thermally-regulated envelope phase (see Section 2.4) shown 
in the left panel of Figure 1 . Such rapid migration rates may 
be difficult to reconcile with formation via core accretion 
of giant planets that orbit at large distances from their stars, 
unless the time spent by a planet in the slow growth stage 
and the migration timescale t\ are actually closer (see, e.g., 
Alibert et al, 2005). 

Once a planet becomes a gas giant, a density gap forms 
along its orbit (because condition 40 is typically satisfied) 
where the surface density is drastically depleted (see right 
panel of Figures 3 and 5). In the limit of a very clean gap, 
there is a balance between viscous torques and gravitational 
torques at gap edges (interior and exterior of the planet's or- 
bit) and the planet remains locked in the gap. This argument 
is valid as long as the disk inside the orbit of the planet is 
not significantly depleted. The planet drifts towards the star 
carried by the viscous diffusion of the disk, if the local mass 
of the disk is comparable to the planet's mass, or larger. 
This is referred to as type II migration. Viscous diffusion 
through the disk occurs on a timescale a \ da/dt \ ~ 1 ~ a 2 / v , 
and therefore 

^„~i(!) 2 . (47) 

In typical circumstances, there is residual gas in the 
gap region and there may not be an exact balance between 
viscous and gravitational torques at gap edges. Nonethe- 
less, results from multi-dimensional simulations indicate 
that migration occurs on a timescale of order given by equa- 
tion (47) in the presence of a sufficiently deep gap (e.g., 
DAngelo and Lubow, 2008). In a disk with a ~ 10 -3 
and H/a ~ 0.05, a gap-opening planet would undergo or- 
bital migration on timescales of many tens of thousands of 
orbital periods, or many 10 5 years at ~ 5 AU from a solar- 
mass star. 

According to equation (47), m oc a 3 / 2 (if H/a is ap- 
proximately constant) and thus the migration timescale be- 
comes shorter as a giant planet approaches the star. How- 
ever, efficient exchange of orbital angular momentum be- 
tween a gas giant and a protoplanetary disk requires that 
the local disk mass is at least comparable with the planet's 
mass, i.e., 7ra 2 E > M. If this is not the case, planet's in- 
ertia acts to slow down migration. Strict type II migration 
of a 1 Mj planet would therefore necessitate a gas surface 
density £ at both gap edges greater than ~ 100 gem -2 at 
5 AU from the star and > 2.8 x 10 3 g cm" 2 at 1 AU. 

If the disk inside the planet's orbit is significantly de- 
pleted (due to the tidal barrier of a massive planet), conser- 
vation of angular momentum at the gap edge exterior of the 
planet's orbit requires that dL p /dt + % = (e.g., Syer and 
Clarke, 1995), where the first term on the lefthand side is 
the rate of change of the planet's orbital angular momentum 
and the second term is the (positive) viscous torque exerted 
on the disk, given by equation (39). Thus, for a nearly- 
truncated disk at the gap edge exterior of the planet's orbit 



da/dt = —2% / (MaQ) , and the migration timescale is 

°» ~ L (f ) 2 {£) ■ <«•> 

Inertia effects are expected to be important at small orbital 
radii, when a disk is sufficiently depleted. Migration of gi- 
ant planets provides a natural explanation for the occurrence 
of Jupiter-mass planets orbiting within a few tenth of AU 
from their parent stars, and whose existence would be oth- 
erwise difficult to explain by means of in situ formation via 
core accretion or disk instability (as we shall see in Sec- 
tion 3). 

2.10. Examples of Giant Planet Evolution Tracks 

An illustrative example of models combining planet's 
growth and orbital migration is displayed in Figure 6. These 
results are obtained from calculations that encapsulate in 
a simplified manner all the basic aspects of a core accre- 
tion model and standard migration theory discussed thus 
far. We consider here the case of a planetary embryo em- 
bedded in an axisymmetric (around the star) disk, which 
evolves under the action of viscous torques (equation (39)), 
gravitational torques, and photoevaporation induced by the 
central star. All three effects are accounted for in equa- 
tion (41) once the lefthand side of that equation is replaced 
with 7rad(H + £ pe ) jdt , in which d£ pe /<% is the mass loss 
flux due to photoevaporation. 

Tidal torques exerted by the planet on the disk are incor- 
porated in equation (41) through a set of torque density dis- 
tributions (torque per unit disk mass as a function of radius), 
which cover the relevant range of planet masses and are ob- 
tained from three-dimensional hydrodynamical simulations 
of disk-planet interactions (see DAngelo andLubow, 2008, 
for details). This procedure allows us to encompass the 
different regimes of orbital migration experienced by the 
planet as it grows and discussed in Section 2.9. 

We use a simplified form of the mass loss rate from 
the disk that is zero within the "gravitational" radius 
a 9 w 10 (M*/M ) AU, where the thermal energy of 
the gas equals in magnitude its gravitational energy, oth- 
erwise dT> pe /dt = T,^ e (a g /a) 5 ^ 2 (see the chapter by 
Roberge and Kamp). The mass loss flux at radius a g is 

E£ e = 3.7 x 10- 13 ($pi) 1/2 (M*/M )~ 3/2 in units of 
Mq AU -2 per year. The constant ^\ is the rate of ioniz- 
ing photons emitted by the star, <£ pe , in units of 10 41 s _1 . 

At time t = 0, the disk contains « 0.022 M of gas 
within about 40 AU (Davis, 2005). The disk evolution 
equation is solved numerically, adopting a turbulence pa- 
rameter ft^4x 10 -3 and a disk scale-height H/a = 0.05. 
The evolution begins at 13 AU from a solar-mass star. 

A planet embryo, of initial mass 0.1 M , accretes solids 
at an oligarchic rate (oc M 2 / 3 , see Section 2.2) and be- 
comes isolated when M « 13 M . The gas accretion rate, 
dM/dt, is given by the smaller of M/r e (equation (15)) 
and the upper fit to data in Figure 4. Three values of 
are used: 1, 10, and 50. Thicker and thinner lines represent 
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Figure 6. Evolution tracks of giant planet formation via core accretion obtained from models that account for orbital migration by tidal 
interactions with the protoplanetary disk (see Section 2.10 for details). The total mass (left) and orbital radius (right) of the planet are 
shown as a function of time. Two sets of tracks represent cases in which the planet evolution begins at time t = (thicker lines) and 
at time t = 1.5 Myr (thinner lines). Within each set of tracks, three different line types represent different values of the applied rate of 
ionizing photons emitted by the star: $ pe = 10 41 s _1 (solid lines), 10 42 s _1 (long-dashed lines), and 5 x 10 42 s _1 (short-dashed lines). 



cases in which the planet's evolution starts at time t = 
and 1.5 Myr, respectively. The evolution tracks illustrated 
in Figure 6 highlight the importance of the disk's initial con- 
ditions in which a giant planet's core begins to evolve and 
on the details of disk photoevaporation. It also supports the 
contention that early core accumulation may lead to sub- 
stantial amounts of radial migration, whereas late- stage for- 
mation scenarios may imply much less orbital decay. 

2.11. Summary of Core Accretion Models 

The formation of a gas giant via core nucleated accre- 
tion may take about one to a few Myr and is initiated by 
a solid planetary core of at least a few Earth masses. Core 
formation requires that the heavy elements of the disk have 
condensed out and coagulated into planetesimals. Although 
the process of planetesimal formation is still under scrutiny, 
small bodies of ~ 10 to ~ 100 km in size are observed 
in the Kuiper belt and Neptune-mass planets are observed 
around the Sun and other stars. Core accretion models pro- 
vide an explanation for the existence not only of the four 
giant planets in the solar system, but also for the majority 
of giant planets that have been observed around other stars. 

Giant planets contain large amounts of Hydrogen and 
Helium, which are acquired in gaseous form from proto- 
planetary disks. Therefore, gas dispersal timescales of disks 
set a strict upper limit on the time available to form these 
planets (a few to ~ 10 Myr). The phase of slow envelope 
contraction that precedes runaway gas accretion seems to 
represent the greatest hurdle to overcome since it may last 
longer than disk's gases, at least in some situations. Orbital 
decay by tidal torques may require formation in evolved 
rather than young disks to avoid amounts of radial migration 
that are too large. Additionally, since the oligarchic growth 
timescale is oc a?/ 2 M* 1 / 2 , core formation at orbital dis- 
tances of many tens of AU or around low-mass stars (e..g., 
red dwarfs) may take too long to produce a gas giant. 



3. DISK INSTABILITY MODELS 

Stars form through the gravitational collapse of inter- 
stellar clouds, and Boss (1997) has championed the idea 
(Kuiper, 1951; Cameron, 1978) that gas giant planets might 
form from protoplanetary disks in a similar way. This is 
commonly referred to as the disk instability theory of gas gi- 
ant planet formation, because the mechanism by which the 
planets form is the onset of gas-phase gravitational instabil- 
ities (hereafter referred to as GIs) leading to fragmentation 
of a protoplanetary disk into bound, self-gravitating clumps. 
An important distinction with core accretion is that, if the 
conditions prevail for disk instability to work, then planets 
can form directly out of the gas phase within only a few 
to tens of disk orbit periods (r ro t = 2ir/Q). The planets 
gain most of their complement of gas immediately and sub- 
sequently sediment heavy element cores and/or sweep up 
planetesimals over a longer timescale. Disk instability is 
top down and initially rapid; core accretion is bottom up 
and initially slow. The critical issue for disk instability is 
when and where the conditions necessary for fragmentation 
might actually occur in protoplanetary disks. 

We will begin by discussing how GIs operate and what 
is required for the outcome to be fragmentation. To decide 
whether real disks fragment, we must rely heavily on de- 
tailed numerical simulations. 

3.1. Gravitational Instabilities: Linear Instability 

Gravitational instabilities are caused by the self-gravity 
of the disk, but, unlike star formation, the onset of planet 
formation via disk instability is much more strongly influ- 
enced by rotation, and it happens in a different thermody- 
namic environment. 

To understand how GIs work, consider first the basic 
equilibrium of a gas disk orbiting a central star for the case 
where disk self-gravity cannot be ignored. For the vertical 
equilibrium in a thin disk where parameters vary smoothly 
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and slowly with orbital distance a, equation (3) becomes 

where E(a, z) = f* pdz' is the surface density of gas at 
radial position a between — z and z, with z = being the 
midplane. In what follows, we use E = E(a) to refer to 
E(a, oo ). The additional (third) term on the lefthand side 
of equation (49) characterizes the vertical component of the 
gravitational field due to the gas disk. The radial equilib- 
rium for a thin disk is then 



1 d P n 2 
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0, 



(50) 



where is the additional radial component of the grav- 
itational field due to disk self-gravity, which depends on 
the detailed distribution of the disk surface density E(a). 
When H/a is small, disk self-gravity can be important in 
equation (49) without necessarily being a major contribu- 
tion in equation (50). So, to derive simple relations in later 
sections, we will sometimes approximate Q by using the 
Keplerian angular speed, y/GM^/a 3 . 

GIs arise through perturbation of the equilibrium quan- 
tified by equations (49) and (50). A classic linear sta- 
bility analysis by Toomre (1981) for axisymmetric radial 
ripples is very instructive. Equations (49) and (50) de- 
fine the zero-order equilibrium state of the disk, which 
we denote by a subscript zero, e.g., po, po, etc. Physi- 
cal variables in the perturbed disk can then be represented 
by these zero-order quantities plus small perturbations, i.e., 
P = Po + tip, with \Sp\ <C po- To simplify the analy- 
sis, the radial wavelength A of the perturbations is taken 
to be small compared with both a and the scale length of 
radial changes in E, but large compared with the vertical 
scale-height H. The thin vertical structure of the disk is 
assumed to remain in hydrostatic equilibrium according to 
equation (49) even when perturbed. So, after the accelera- 
tion term dv a /dt ~ d(Sv a )/dt is introduced into the left- 
hand side of equation (50) in order to describe the dynamics 
of the perturbations, the radial equation is integrated over z. 
Equation (50) allows the zero-order terms to be cancelled, 
and only terms to first order in the perturbed quantities are 
retained in the linear analysis (for more details, see Binney 
and Tremaine, 1987). 

By looking for wave-like perturbations such that, e.g., 
Sp oc exp (iut ± 2ira/\) , one can derive the dispersion re- 
lation 

2 2 4^ 2 GE 4^ 2 c 2 
u 2 = n 2 — + — — , (51) 

where A is positive-definite and c 2 is the sound speed 
squared. The quantity k in equation (5 1) is the epicyclic fre- 
quency at which a perturbed fluid element oscillates about 
its equilibrium position {Binney and Tremaine, 1987). For 
a strictly Keplerian disk, k, = Q, because, when the radial 
equilibrium is dominated by the central force of the star, a 
slightly perturbed fluid element follows an elliptical orbit 



with the same period as r ro t- As viewed from its unper- 
turbed equilibrium position, an oscillating fluid element ap- 
pears to execute a small elliptical "epicycle". The epicycle 
has the opposite or retrograde sense of rotation compared 
with the orbital motion of the disk. 

Equation (51) reveals a great deal about the behavior of 
perturbations in self-gravitating disks. As long as the right- 
hand side (RHS) of equation (51) is positive, iut ± 27ra/A 
is a complex number and so the ripples are waves which os- 
cillate and propagate radially. However, when the RHS is 
negative, there are stationary solutions that grow exponen- 
tially (the argument of the exponential is real). The disk is 
then unstable and tends to break apart into dense concentric 
rings. 

An examination of the individual terms on the RHS of 
equation (51) reveals the stabilizing and destabilizing in- 
fluences. The first and third terms are always positive and 
represent stabilization by rotation and gas pressure, respec- 
tively. As A gets large, the first term dominates the RHS of 
equation (5 1); in other words, rotation stabilizes long wave- 
lengths. Because the central gravitational force of the star 
appears in equation (51) through k, it is also a stabilizing 
influence. When A is small, the third term dominates; in 
other words, gas pressure stabilizes small wavelengths. In 
contrast to the first and third terms, the self-gravity of the 
disk represented by the second term on the RHS , which is 
only important at intermediate wavelengths, is always neg- 
ative and hence is a destabilizing influence. If the accel- 
eration imparted by disk gas, G E, is large enough relative 
to ft 2 A and c 2 /A, some middle range of wavelengths will 
be unstable. Because the destabilizing term is due to disk 
self-gravity, we refer to such instabilities as gravitational 
instabilities. 

Setting uj = in equation (51) and solving for A, one 
can see by simple algebra that there is a range of unstable 
wavelengths when 



ttGE 



< 1. 



(52) 



Notice that the numerator of the Toomre stability param- 
eter Q is, in effect, a geometric mean of the stabilizing in- 
fluences of pressure and rotation, while the denominator 
represents the destabilizing influence of disk self-gravity. 
Thus, fixing the mass of the star, a protoplanetary gas disk 
becomes unstable to GIs if it is sufficiently cold (low c) or 
massive (high E). The critical wavelength A cr that first be- 
comes unstable when Q dips below unity is 



Apr — 



2tt 2 GE _ 2ttH 



(53) 



as can be derived by setting w = and Q = 1 in equa- 
tion (51). For the approximate relation on the right in equa- 
tion (53), we use c w HQ, which still holds to within a 
factor of order unity for a self-gravitating disk. 

As summarized in various reviews (e.g., Durisen et aL, 
2007), multi-dimensional simulations of disks show that 
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nonaxisymmetric perturbations, i.e., spiral waves, become 
unstable for 

Q=^<Qcr, (54) 

where Q cr is between about 1.5 and 2.0, depending on the 
detailed structure of the gas disk. Because spiral waves can 
be unstable in disks that are stable against growth of the 
axially symmetric ripples we just discussed, spiral distur- 
bances are considered to be more unstable, i.e., a disk which 
transitions slowly from stable to unstable conditions will 
first manifest GIs in the form of nonaxisymmetric waves. 
In fact, in simulations, the linear growth of GIs from noise 
is characterized by the appearance of multi-armed spirals, 
which are typically trailing due to the Keplerian shear. 
These waves grow exponentially from noise in the linear 
regime on the dynamic timescale r ro t (e.g., Nelson et al., 
1998; Pickett et al, 1998). 

With a little algebra, Q in equation (52) can be rewrit- 
ten approximately as twice the ratio of the star's vertical 
gravitational field to the disk's vertical gravitational field, 
namely, 2g*/g® , where both g\ and are evaluated at 
z = H. So, Q < 2 corresponds roughly to g^ > g*. 
As expected, GIs appear when the self-gravity of the gas 
disk becomes sufficiently important, i.e., when the second 
and third term on the lefthand side of equation (49) become 
comparable. 

The discussion associated with equation (51) makes it 
clear how self-gravity could break a disk into ringlets, but 
how do spiral waves grow? The likely mechanism is ex- 
plained in Toomre (1981) and Binney and Tremaine (1987). 
Consider part of a leading spiral wave spanning its own 
corotation radius (CR hereafter), the radial location in the 
disk where the angular pattern speed of the wave equals ft, 
i.e., where one arm of the wave goes around the disk in ex- 
actly r ro t • Shear in the disk causes the arc to "swing" from 
leading to trailing. Because k ~ ft in a nearly Keplerian 
disk, the fluid elements of the disk in the crest of the wave 
near the CR are executing epicyclic oscillations that keep 
them in the crest. This allows self-gravity to amplify the 
density maximum at the crest of the wave. Numerical stud- 
ies of protoplanetary disks in which spiral waves grow show 
that the conditions are indeed right for the swing amplifi- 
cation mechanism to operate (Pickett et al., 1998; Mayer 
et al., 2004), and the spiral waves that grow always straddle 
their CR (see. e.g., Mejia et al., 2005). 

3.2. Occurrence of Gravitational Instabilities 

To get a feeling for when GIs are expected in protoplan- 
etary disks, we can evaluate equation (54) for the typical 
parameters used in earlier sections, namely H/a w 0.05 
at a = 5.2 AU for a disk orbiting a solar-mass star (corre- 
sponding to a temperature T ~ 75 K if c is the adiabatic 
sound speed). Let us assume that Q cr ^1.7. Then, 

q*q„( 3xl °; 8Cm - 2 ). (55, 



The typical solids' surface density E s ~ 10 gem -2 dis- 
cussed in Section 2 for core accretion models scales to 
E ~ 10 3 g cm -2 for gas with solar composition when half 
the heavy elements are assumed to be in the solid phase. 
A surface density £ = 10 3 g cm -2 at 5.2 AU is three 
times the standard Hayashi MMSN (Hayashi, 1981), but 
is only a factor of about 1.5 larger than more recent MMSN 
models (Davis, 2005; Desch, 2007). For a disk around a 
solar-type star to be unstable to gas-phase GIs at Jupiter's 
orbit radius, i.e., Q < Q cr , equation (55) tells us that E 
must be > 3 x 10 3 g cm -2 , which is ten times higher than 
the Hayashi MMSN and four or five times higher than the 
newer MMSN models. 

There are two ways that such a high E might occur at 
5.2 AU: either the disk is massive everywhere or there is 
a localized concentration of mass near this radius. In the 
first case, where E(a) is relatively smooth, gas-phase GIs 
at the orbit radii of the gas giants in the solar system re- 
quire a total protoplanetary gas disk mass of about 0.1 M 
or more (Boss, 2002a), and disk masses of this order are 
typical in hydrodynamics simulations used to study GIs in 
protoplanetary disks. Disks may well have masses this large 
or larger during the accretion phase, when gas is falling onto 
the young star/disk system at rates of 10 -5 M per year or 
higher (e.g., Vorobyov and Basu, 2006). Alternatively, a 
disk might have a localized enhancement of surface den- 
sity due to accumulation of mass in or at the edges of a 
dead zone (Gammie, 1996), where radial transport of mass 
becomes inefficient. This could lead to episodic eruptive 
phenomena related to FU Orionis outbursts (e.g., Armitage 
et al., 2001; Boley and Durisen, 2008; Zhu et al., 2009). 

Although we have not yet discussed the conditions for 
fragmentation, it is instructive to estimate how much mass 
is involved in the GIs. Using equation (53) as an estimate 
for the radial extent of the instability, we expect the mass in 
one arm of a /c-armed unstable wave at 5.2 AU to be about 

M frag «2 7 r^~^M J . (56) 

Simulations show that k is typically four or five at the onset 
of GIs (Mejia et al., 2005; Boley and Durisen, 2008) when 
Mjj ~ 0. 14 M* , and therefore the arms of the spirals would 
indeed have a gas inventory sufficient to produce a giant 
planet. 

Even if a disk is stable (Q > Q cr ) at the radius of 
Jupiter's orbit, it may be susceptible to GIs at larger radii. 
Suppose that the disk temperature T cx a -1 / 2 (i.e., H/a cx 
a 1 / 4 , a "flared" disk). We know that k « ft cx a -3 / 2 , so if 
E cx a~ s , 

Thus, for any disk with s < 7/4, Q becomes less than Q cr 
if the disk extends to a sufficiently large radius. For a disk 
with s < 2, M D (a) = /^WEda' ~ a 2 ~ s , and so the 
disk mass Md(cl) interior of a diverges as a increases. As a 
result, if a disk with s < 7/4 is extended enough to have GIs 
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in its outer regions, it is likely to have a mass that is a signif- 
icant fraction of the star's mass. Again, this is most likely to 
occur during the early accretion phase of the star/disk sys- 
tem. If we put in parameters appropriate for this situation 
at a few hundred AU around a solar-type star, Mf rag can 
become tens of Jupiter's masses or more (Stamatellos and 
Whitworth, 2009). 

The above arguments are not very sensitive to the tem- 
perature distribution T(a). If we use T(a) ex a or oc a -1 
instead, then we get H/a oc a 1 / 2 or H/a equal to a con- 
stant, respectively, and Q oc a s_3 / 2 or Q oc a s ~ 2 , respec- 
tively. Parameter Q will still become less than Q CY at large 
enough radii for surface densities that do not fall off too 
rapidly. During the early accretion phase, one expects that, 
in spatially extended disks, the gas temperature levels off 
to a constant value (T oc a ) at large radii due to envelope 
irradiation. At small radii, some disks may not be irradiated 
by their central stars due to shadowing. Such disks would 
have steeper temperature f ail-off s than a -1 / 2 and may not 
be flared. For more detailed analytic arguments about grav- 
itational instabilities in disks during the accretion phase, see 
Clarke (2009) and Rafikov (2009). 

3.3. Nonlinear Growth of Gravitational Instabilities 
and Fragmentation 

There are two possible nonlinear outcomes for GIs: 
steady- state balance of heating and cooling or fragmenta- 
tion into dense clumps. As the spiral waves grow in ampli- 
tude and steepen into shocks, they turn some of the ordered 
rotational energy of the disk into heat. Moreover, gravi- 
tational torques due to the predominantly trailing charac- 
ter of the spirals lead to net outward transport of angular 
momentum and inward transport of mass. This transport 
causes gravity to do work on fluid elements as they sink 
into the central gravitational potential and is another impor- 
tant source of heat for the gas. If the disk does not shed this 
heat, its temperature and hence c inevitably becomes larger, 
which raises Q and shuts off the instability. Disks do radi- 
ate, and so GIs can be sustained in an quasisteady state with 
an average Q near but somewhat below Q cr where the heat- 
ing by GIs is balanced by radiative cooling (Gammie, 2001 ; 
Lodato and Rice, 2004; Mejia et al., 2005). Mejia et al. re- 
fer to this as the asymptotic state. GIs grow on a dynamical 
timescale r rot for Q well below Q cr , but the growth times 
become long as Q — > Q cr from the unstable side. So the 
heating rate due to GIs can adjust to the prevailing radiative 
cooling rate. In this sense, the thermal physics of the disk 
controls the limiting amplitude of GIs (Tomley et al., 1991; 
Pickett et al., 2000). 

Let us define a radiative cooling time in a column- wise 
sense by 



7~cool 



(58) 



2 0~ rad^eff 

where E[ nt is the gas internal energy density and the effec- 
tive temperature, T e ff , characterizes the radiant flux (energy 
per unit area per unit time) a ra dT e 4 ff out the top of a col- 



umn through the disk. The factor of two in equation (58) 
accounts for radiation from both the top and bottom of the 
disk. This equation makes the approximation that all en- 
ergy loss from the disk is vertical, which need not be true in 
complicated dynamical situations. 

When t coo i is comparable to or less than r ro t , the balance 
of heating and cooling discussed above may not be possi- 
ble. The gas could cool profoundly between the spiral arms 
before it is reheated by the shock in the next arm. Another 
way to say the same thing is that, because the growth time of 
GIs is r ro t , it is the shortest timescale on which GIs can re- 
plenish thermal energy in the gas. There will then be large 
compressions behind the shocks where the self-gravity of 
the gas could win out entirely. The arms may then break 
up into dense clumps, which become bound in the sense 
that the magnitude of the clump self-gravitational energy 
exceeds the clump's own internal and rotational energies. 

The precise value of r coo i for which dense clumps form 
must be determined numerically. The first systematic treat- 
ment by Gammie (2001) for a simple equation of state in a 
local patch of a razor thin disk showed that fragmentation 
occurs for ^r coo i < 3, or equivalently r coo i/r rot < 1/2. 
Global three-dimensional simulations by several groups us- 
ing a variety of simulation techniques have verified that the 
fragmentation criterion is 



7~cool 
7"rot 



< 



(59) 



where /f rag depends on the gas equation of state and, if 
T coo i is a function of time, on the thermal history of the 
gas (Rice et al., 2003, 2005; Mejia et al., 2005; Clarke 
et al., 2007). For ideal gases with ratios of specific heats 
T\ = 5/3 and 7/5, /f rag « 1 and 2, respectively. In other 
words, a "softer" equation of state in the sense of lower Ti 
makes a disk more unstable. This may be relevant, because 
5.2 AU is close to where the effective Ti for molecular Hy- 
drogen is expected to change from 5/3 at large disk radii 
to 7/5 at intermediate radii. When r coo i is time-dependent 
and approaches the fragmentation condition slowly, /f rag 
becomes about half as large as when r coo i is constant, and 
so the disks become more stable and harder to fragment. 

For the special case of gas disks that behave isother- 
mally when perturbed, r coo i is essentially zero, and, instead 
of a fragmentation criterion based on the ratio r coo \/r YOt , 
one finds that all disks fragment when Q < Qf rag where 
Qfrag ~ 1.4-1.5 < Q cr (Boss, 2000; Mayer et al., 2002; 
Johnson and Gammie, 2003). A classic example of frag- 
mentation in an isothermally evolved disk is shown in Fig- 
ure 7. 

The masses of planets that form in fragmenting simu- 
lations are roughly compatible with equation (56), but the 
process of fragmentation can be violent and chaotic, of- 
ten leading to multiple interacting protoplanetary clumps. 
Some clumps merge or shear out, and new clumps can form. 
No one has yet integrated a large number of fragmenting 
simulations over long times, but a moderately broad spec- 
trum of final clump masses would probably result in a mass 
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Figure 7. Midplane density contours after 374 years for an 
isothermal simulation of a protoplanetary gas disk spanning 4 to 
20 AU around a solar-mass star with Md — 0.09 M and an 
initial Q m in = 1.3. The dense clump near twelve o'clock is grav- 
itationally bound and has about 5 Mj . Adapted from Boss (2000). 

range between tenths to multiples of Mj for disks with 
Md ~ 0.1 M* that are a few tens of AU in radius (see, 
e.g., Mayer et al., 2004). A discussion of fragmentation 
spectra for disks with more extended radii can be found in 
Stamatellos and Whitworth (2009). 

3 .4 . Realistic Radiative Cooling 

The pertinent question now is whether r coo i for realistic 
dust opacities is actually short enough to cause fragmenta- 
tion. The vertical optical depth down into a disk is defined 
by 

r(z) = / Rpdz' , (60) 

J z 

where R, with units of cross-section per unit mass, is the 
mass absorption coefficient for radiation (opacity) averaged 
in an appropriate way over frequency. If R does not vary 
much with z, then the midplane optical depth is r(0) ~ 
RTt/2. This is essentially the number of photon mean free 
paths between the midplane and z = oo . 

If the radiation emitted by the gas and dust is thermal and 
if a frequency-independent (grey) opacity characterized by 
R is a good approximation, then momentum conservation 
for radiation gives (e.g., Gray, 1992) 

%^ = ( T^) F«* ^ 

dr V 4cr rad/ 

for plane-parallel geometry, where p ra d is the radiation 
pressure, F ra d is the net upward flux of radiant energy, 
and a ra d — 7.57 x 10 -15 erg cm -3 K -4 is the radiation- 
density constant (notice that the ratio in parentheses in 
equation (61) is equal to the inverse of the speed of light). 

For no radiation shining down on the disk, for a con- 
stant upward flux (F ra d = cr ra dT e 4 ff ), and for a radiation 
field that is well described by black-body properties, e.g., 



Pr&d = & m d^ 4 /3, the standard grey atmosphere solution to 
equation (61) is 

T 4 = 3 -T^(r + q), (62) 

where q depends on boundary conditions and approxima- 
tions but is roughly 2/3. Equation (62) strictly applies only 
to a semi-infinite atmosphere with a constant upward flux. 
A disk is finite and has a flux which varies with r (Hubeny, 
1990). Nevertheless, within a factor of order unity, equa- 
tion (62) is a useful estimate relating the midplane temper- 
ature T mid = T(z = 0), i.e., T(r) at r = r(0), to the 
expected effective temperature T e ^ in radiative equilibrium 
given a midplane optical depth r (0) . 

When r(0) ^> q, the disk is optically thick. In equa- 
tion (58), one can approximate the integral of E{ nt to be 
3£7£T m id/(2/i), where 1Z is the gas constant and /i is the 
mean-molecular weight (see Section 2.5). Then, 

where, for simplicity, the final dependence on a is based 
on assuming that R is constant and T m id oc a -1 / 2 . (If 
T m id oc a or oc a _1 then r coo i oc a~ 2s or r coo i oc a 3_2s , 
respectively). At 5.2 AU for the usual parameters, equa- 
tions (60) and (63) give 

and 

r cool ~10 5 yr(- ^—j.) (— -| -s) . (65) 

\lcm^g l J \3xl0 d gcm 1 ) 

Note that, for H/a = 0.05 at 5.2 AU, T mid « 75 K if we 
use an adiabatic sound speed with I\ = 5/3 and /i w 2.4. 
In an unstable disk at 5.2 AU, where r rot ~ 12 years, equa- 
tion (65) gives T C ooi/Vrot ~ 10 4 , which is orders of magni- 
tude greater than /f rag for any reasonable value of Ti . 

We have just seen that, for our standard parameters, 
equations (54), (55), (59), and (63) together imply that, at 
5.2 AU, a disk which is dense enough to be gravitationally 
unstable will cool too slowly to fragment. Rafikov (2005) 
first pointed out this conundrum. Although it depends on 
detailed assumptions about how disk properties vary with a, 
it can reasonably be argued that this problem extends out to 
about 40 AU or so, exactly the region where one would like 
to form the giant planets in our own solar system. More- 
over, if we relax our choice of H/a = 0.05, then lowering 
T coo i at fixed R requires a lower E and/or a higher T m id, 
both of which would tend to make the disk more stable. 

One can instead imagine lowering R itself because dust 
grows and settles to the midplane, but there is a limit to the 
efficacy of this ploy. The most efficient cooling by thermal 
emission occurs for r(0) ~ 1. For lower r(0), the disk 
is a poor absorber and, hence, by Kirchhoff 's law, a poor 



D'Angelo, Durisen, & Lissauer 



19 of 28 



Giant Planet Formation 



3 .5 Computational Results for Disk Instability 



3 DISK INSTABILITY MODELS 




Figure 8. Midplane density grey scale of a nonfragmenting and 
radiatively cooled disk with M* = 0.5 M , M D = 0.07 M , 
and an initial Q m in = 1.5 after about 2000 years. The box shown 
is 80 AU on a side. Adapted from Boley et al. (2006). 

thermal emitter. All else being equal, equation (64) shows 
that a reduction of r(0) to unity requires a reduction of R 
by a factor of 10 3 . With this reduction, r coo \/r rot becomes 
about 10, which is still larger than the largest expected value 

Of /frag- 

These analytic arguments suggest that gas giant planet 
formation by disk instability probably does not work in- 
side several tens of AU. On the other hand, equations (57) 
and (63) further suggest that the instability criterion (equa- 
tion (54)) and the fragmentation criterion (equation (59)) 
may be satisfied simultaneously at large enough radii, 
roughly 100 AU and beyond, when < s < 7/4 (see 
also Clarke, 2009; Rafikov, 2009). 

3.5. Computational Results for Disk Instability 

3.5.1 . Simulations with Realistic Radiative Cooling 

The assumptions that go into the results of Section 3.4 
are rather simplistic and the details are subject to debate. 
For instance, R is a complicated function of particle size and 
temperature and T m id at very large orbital radii may reach 
a constant value. It can also be tricky to define r coo i prop- 
erly for the fragmentation criterion in a radiatively cooled 
disk evolving under the action of GIs {Johnson and Gam- 
mie, 2003). Disk instability must be tested by detailed mul- 
tidimensional radiative hydrodynamics simulations. The 
numerical techniques for such simulations remain difficult, 
and results to date are disparate and controversial. A de- 
tailed discussion of numerical schemes goes beyond the 
scope of this chapter. Here we only summarize the current 
situation regarding global simulations that attempt to treat 
the cooling correctly with realistic opacities. This section 
deals with simulations of disk instability in regions extend- 
ing out to tens of AU. 

The analytic arguments of Section 3 .4 suggest that frag- 



mentation should not occur inside 40 AU or so. Evi- 
dence supporting this conclusion comes from simulations 
using both Smoothed Particle Hydrodynamics (SPH) (Nel- 
son et al., 2000; Stamatellos and Whitworth, 2008; For- 
gan et al., 2009) and grid-based hydrodynamics (Cai et al., 
2006; Boley et al, 2006, 2007; Boley and Durisen, 2008). 
These papers employ a variety of radiative algorithms. Fig- 
ure 8 illustrates a typical midplane density structure for a 
nonfragmenting simulation. These studies treat 0.5-1 M 
stars with a disk mass Mjj ~ 0.1 M* and with outer disk 
radii of ten to tens of AU. The disks do not cool rapidly 
enough to fragment. In the cases where r coo i is explicitly 
calculated, nonfragmentation is consistent with the appli- 
cation of equation (59) for the relevant equation of state. 
Some efforts have been made by those who do not see frag- 
mentation to evolve similar disks and to test their radiative 
schemes against analytic expectations, e.g., compare Boley 
et al. (2007) and Stamatellos and Whitworth (2008). 

Evidence for fragmentation comes from two groups, one 
using grid-based simulations (Boss, 2001, 2002a,b, 2007, 
2008) and the other using SPH (Mayer et al., 2007). The 
disk parameters are similar to those of the nonfragmenting 
simulations, but these researchers find that their disks frag- 
ment into dense, bound clumps, with masses of typically a 
few Mj (see equation (56)). These groups do not agree in 
detail on the conditions necessary for fragmentation. Boss 
finds that fragmentation is robust, i.e., fairly independent 
of many physical and numerical parameters, while Mayer 
et al. find that fragmentation requires somewhat higher disk 
masses and high mean molecular weight and that it is sen- 
sitive to some choices of numerical parameters. Cai et al. 
(2010), using an improved radiative cooling scheme simi- 
lar to that of Boley et al. (2006) but with higher numeri- 
cal resolution, have now evolved the same initial disk and 
perturbation as in one of the Boss (2007) simulations. No 
fragmentation occurs. The physical conditions are almost 
identical, so Cai et al. (2010) conclude that the difference 
really seems to be a matter of radiative algorithms. 

The different outcomes are not due to inadequate resolu- 
tion in the nonfragmenting calculations. For instance, some 
of the simulations in Boley and Durisen (2008) have ex- 
tremely large grids, achieving resolutions of 0.025 AU over 
a radial range of about 13 AU. This resolution is more than 
adequate to detect fragmentation when conditions are arti- 
ficially forced to satisfy equation (59). 

Both Boss (2004) and Mayer et al. (2007) recognize that 
purely radiative cooling is not fast enough to satisfy equa- 
tion (59). Both attribute the fragmentation they see to rapid 
cooling due to convective transport in the vertical direction 
and cite upwellings at the sound speed associated with spi- 
ral arms in their disks as evidence for thermal convection. 
However, numerical studies (Boley et al., 2006, 2007) com- 
bined with analytic arguments (Rafikov, 2007) suggest that 
the grid-based codes which do not see fragmentation are 
fully capable of modeling thermal convection and that ther- 
mal convection in a disk is unlikely to carry more than tens 
of percent of the vertical flux, insufficient to reduce cooling 
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times by orders of magnitude. Boley and Durisen (2006) 
explain the upwellings associated with spiral shocks, which 
are also seen in nonfragmenting simulations, as hydraulic 
jumps, a dynamic phenomenon very different from thermal 
convection. 

The disagreement about cooling times leads to disagree- 
ment about how the disk metallicity affects fragmentation. 
Boss (2002b) studies variations in metallicity by varying R 
up and down by factors of ten relative to solar composition 
and finds that fragmentation is not sensitive to R. On the 
other hand, Cai et al. (2006), using simulation techniques 
similar to those in Boley et al. (2006), find that GIs are sen- 
sitive to metallicity in simulations where R is varied only 
from 1/4 to 2 times the solar opacity. No simulation in 
Cai et al. produces fragmentation, but the GIs in the lower 
metallicity disks are stronger, which is expected if r coo i is 
mostly determined by radiative processes. 

Finally, Boss (2002a), Cai et al. (2008), and Stamatellos 
and Whitworth (2008) do agree with analytic arguments by 
Matzner and Levin (2005) that sufficiently strong external 
irradiation of a disk can weaken and suppress GIs. Strongly 
irradiated disks should not fragment. An easy way to es- 
timate the necessary flux of irradiation Fi rr = cr ra dX^ r to 
stabilize a disk is to compute the irradiation temperature, 
Ti rr , such that the resulting sound speed c(T- 1TT ) satisfies the 
inequality in equation (54). At 5.2 AU, T- 1TT ~ 75 K suffices 
to suppress GIs altogether for Mb ~ 0.1 M . No matter 
what else one might believe about disk instability in the re- 
gion within tens of AU, it seems pretty clear that irradiation 
of this order or higher prevents planets from forming by 
GIs. 

3.5.2. Enrichment of Heavy Elements 

In the standard disk instability model, heavy element 
cores form after the bulk of the gas in the planet becomes 
a gravitationally bound equilibrium object. Although such 
protoplanets are likely to be rotating (see, e.g., Mayer et al., 
2004), let us approximate them here as hydrostatic equi- 
librium spheres, with the gas pressure gradient balancing 
gravity. As discussed in Section 2.1, solid particles re- 
spond to gas pressure forces only indirectly through gas 
drag. Thus, in the absence of turbulence, small solid par- 
ticles slowly sink to the center of a gas giant protoplanet at 
a drift speed Vd ~ grf (see Section 2.1), where g is the 
magnitude of the gravitational acceleration inside the pro- 
toplanet and Tf is the friction time. For a newly formed 
1 Mj protoplanet with a radius R « 0.5 AU, approxi- 
mately half A cr in equation (53) for our standard parame- 
ters, g « GMj/R 2 ~ 2.3 x 10" 3 cms -2 , and the internal 
temperature would be T ~ 100 K. Using equation (2), one 
gets a sedimentation time 

R _o /lgcm _3 \ / 1 cm \ 
Tsed « 6 x 10 3 yr -2 -5— , (66) 

in rough agreement with an estimate of Boss (1997). 
Detailed calculations by Helled et al. (2008), starting 



with micron- sized particles but including aggregation and 
growth, give similar times. The sedimentation timescale, 
r se d , is short compared with the Kelvin-Helmholtz contrac- 
tion timescale, tkh , of the protoplanet, which is a few times 
10 5 years (Helled et al., 2006). Turbulence slows sedimen- 
tation, but r se d for centimeter- sized particles is still shorter 
than tkh- Therefore, gas giants formed by disk instability 
probably are able to produce heavy element cores out of 
their initial complements of heavy elements. A 1 Mj mass 
of solar composition could form a ~ 6 M§ core. 

Jupiter and Saturn are strongly enhanced in heavy ele- 
ments compared with solar composition, as are some tran- 
siting exoplanets (e.g., HD 149026 b). There are two ways 
that gas giants formed by disk instability could end up 
with complements of heavy elements which are fractionally 
larger than in their natal gas disk: accreting planetesimals 
from the surrounding disk after formation or forming in a 
gas clump already enhanced in solids above the original gas 
disk composition. 

Detailed computations of planetesimal accretion (Helled 
et al., 2006; Helled and Schubert, 2009) for a 1 Mj proto- 
planet formed by disk instability show that total masses of 
accreted planetesimals can vary anywhere from 1 to over 
100 Me, depending upon the initial location of the planet 
and upon the assumed mass and radial density distribution 
of the parental disk. Large accretion rates of heavy elements 
are made possible by the appreciable starting radii of gas 
giants in the disk instability model coupled with their mod- 
erately long contraction timescale. Such calculations are of 
course sensitive to assumptions, and a great deal more work 
is needed, but it appears that heavy element enrichment of 
planets may not be a strong discriminator between disk in- 
stability and core accretion models. 

The other intriguing possibility is that, at the time of 
fragmentation, spiral arms may already be enhanced in 
solids relative to the ambient disk due to gas drag effects. 
How might this occur? We have seen that solid particles 
respond to pressure gradient forces by drifting relative to 
the gas in a direction opposite to the gradient. For exam- 
ple, particles drift or sediment toward the pressure maxi- 
mum at the disk midplane (Section 2.1) or at the center of 
a protoplanet (as discussed above). The same thing applies 
to radial drift of solid particles in the protoplanetary disk: 
they drift radially in a direction opposite to the direction of 
the radial pressure gradient (Weidenschilling , 1977), which 
affect the gas orbital velocity, and concentrate at a radial 
pressure maximum (Haghighipour and Boss, 2003). 

The physics in this case is complicated by angular mo- 
mentum exchange. Equation (50) determines the equilib- 
rium rotation rate of the gas, ft, which includes the effect of 
the gas pressure term. Centrifugal balance against gravity 
in the absence of gas pressure requires a different angular 
speed for solids, ft s , such that ft 2 = GM*/a 3 — g® /a. 
From equation (50) we have that if dp/ da < 0, then 
ft < ft s , while if dp/ da > 0, then ft > ft s . As a 
consequence, gas drag exerts a torque (per unit mass) on 
solid particles equal to (ft — ft s ) a 2 /rf. Then, by conser- 
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vation of angular momentum, the radial drift velocity is 
da/dt = (a/ir) (ft — ft s ) (r rot /rf) (see Section 2.9). In 
disk regions with a negative radial pressure gradient, an or- 
biting particle feels a head wind and drifts inward due to 
loss of orbital angular momentum. In disk regions with a 
positive radial pressure gradient, a solid particle feels a tail 
wind and drifts outward due to gain of orbital angular mo- 
mentum. If there is a pressure maximum at some value of a, 
solids drift toward the pressure maximum from both sides 
and become concentrated relative to the gas. 

The magnitude of the radial drift speed depends on the 
ratio Tf/r rot . For small particles, r//r ro t <C 1 (see equa- 
tion (4)), and particles are essentially entrained with the gas. 
The head and tail winds are very small (ft s c± ft), and so 
the radial drift is slow. Drift speeds increase with parti- 
cle size in the small particle regime. Very large particles, 
for which r//r ro t ^> 1, orbit at angular speed ft s ^ ft. 
They can experience considerable head and tail winds, but 
radial drift speeds are again small because of the long fric- 
tion times. Drift speeds decrease with particle size in this 
regime. Maximum radial drift speeds (~ a\ft — ft s \) occur 
for Tf I r ro t ~ 1 , corresponding to particles of about a me- 
ter in radius (Weidenschilling , 1977). Such particles can be 
swept into radial pressure maxima on timescales of only a 
few hundred years (Haghighipour and Boss, 2003). 

It is not obvious a priori that this phenomenon general- 
izes to the dynamic environment of a Gl-active disk with 
spiral waves, but, in fact, solid particles of this optimal size 
do drift into the pressure maxima of spiral arms on simi- 
larly short timescales (Rice et al., 2003). This concentra- 
tion of solids relative to gas in spiral arms can trigger GIs in 
the solids themselves (Rice et al., 2006) and could even in- 
crease the tendency of the gaseous arms to fragment (Mayer 
etal., 2007). In the latter case, a protoplanet formed by disk 
instability could initially be enhanced in heavy elements rel- 
ative to the disk gas. 

3.6. Bimodal and Hybrid Planet Formation Models 

As argued analytically in Section 3.4 and illustrated in 
Figure 9, even codes that do not produce fragmentation 
in disks within tens of AU around solar-type stars con- 
firm that fragmentation can and does occur in radially ex- 
tended and massive disks, beyond about 100 AU (Boley, 
2009; Stamatellos and Whitworth, 2009). Appropriate scal- 
ing of equation (56) to large orbital distance and very mas- 
sive disks leads us to expect the pieces to be super- Jupiters, 
brown dwarfs, or low-mass stars. This is exciting in light of 
recent discoveries, discussed elsewhere in this volume (see 
the chapter by Traub and Oppenheimer), that massive gas 
giant planets seem to exist at large radii around relatively 
young stars. Standard core accretion models do not form 
planets at large orbital distances because it may take too 
long to form a solid core, whereas violent GIs in an outer 
disk may occur during the early disk accretion phase for a 
range of plausible parameters (Vorobyov and Basu, 2006; 
Clarke, 2009; Rafikov, 2009). 
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Figure 9. Logarithmic surface density (in cgs units) for a radia- 
tively cooled simulation of a fragmenting protoplanetary gas disk 
around a solar-mass star that is accreting at a rate of 10 -4 Mq 
per year. The dense clumps have masses ranging from 4 Mj to 
14 Mj . Adapted from Boley (2009) 

Unless planets formed by core accretion can migrate or 
be scattered to large distances, it seems reasonable to con- 
jecture that planets found in large orbits are formed in situ 
by disk instability. But, if nonfragmentation of Gl-active 
disks in the inner tens of AU prevails as the consensus view, 
we still have to rely on core accretion to explain the bulk 
of gas giant exoplanets discovered so far, whose orbits are 
within a few AU of their stars. This leads to the notion 
(Boley, 2009) that gas giant planet formation may actually 
be bimodal, with different mechanisms dominating in dif- 
ferent regions of protoplanetary disks. Even in this case, it 
is possible that gas-phase GIs play a role in core accretion 
by accelerating the formation of solid planetesimals through 
migration of solids into structures formed in the disk by GIs 
(e.g., Haghighipour and Boss , 2003 ; Rice etal., 2004 , 2006 ; 
Durisen et al., 2005). 

3.7. Summary of Disk Instability Models 

Once the necessary conditions exist, namely Q < 1.5- 
1.7 and T CO oi/T ro t < /frag ^ 1, disk fragmentation into 
bound clumps of gas is a robust dynamic process that re- 
quires only a few to tens of orbital periods. Gas giants 
formed by disk instability can subsequently produce heavy 
element cores by sedimentation and can become enriched 
in heavy elements relative to the background gas disk by 
accreting planetesimals. Disk instability is expected to oc- 
cur in the outer regions of extended massive disks during 
the accretion phase, and it may be the only way to under- 
stand the existence of high planetary-mass objects at large 
distances from their stars, such as Fomalhaut b (Kalas etal., 
2008) and HR 8799 b, c, and d (Marois et al., 2008). 

Although there is still some controversy, analytic argu- 
ments and a number of simulations indicate that disks are 
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unlikely to fragment within a few tens of AU, because the 
conditions of low Q and a low r coo i / T ro t cannot be simulta- 
neously satisfied. The inclusion of stellar irradiation, which 
tends to stabilize these regions by maintaining a high Q, 
makes the problem worse. Disk instability is most likely to 
occur in the outer regions of disks relatively early in disk 
evolution, when disks are massive and still accreting from 
their protostellar cloud at a high rate. The ultimate fate and 
survivability of planets formed at such an early phase of 
disk evolution is unclear. The principle of parsimonious ex- 
planations makes appeal to more than one formation mech- 
anism for gas giants unappealing. 

4. HIGHLIGHTS AND OUTSTANDING QUESTIONS 
4.1. Solar and Extrasolar Planets 

Any theoretical model of giant planet formation has to 
provide a quantitative explanation for the existence of the 
four outer planets of the solar system. Although Jupiter's 
core mass remains uncertain, there is theoretical evidence 
that Saturn, Uranus, and Neptune have cores of 10-20 M e 
(Podolak et al., 1995; Saumon and Guillot, 2004). All four 
planets have massive Hydrogen/Helium envelopes, which 
represent most of the mass of Jupiter and Saturn and more 
than 10% of the mass of Uranus and Neptune (e.g., Guillot, 
2005). 

Core accretion models offer a natural explanation for 
both the formation of heavy element cores and the accre- 
tion of massive gaseous envelopes, such as those of Jupiter 
and Saturn. In this scenario, Uranus and Neptune reached 
their current masses in an environment that may have been 
deprived of gases because their cores took too long to form, 
and never underwent a runaway gas accretion phase. It 
seems unlikely, instead, that any of these planets formed 
via GIs in the early solar nebula, principally because theo- 
retical arguments show that fragmentation may only occur 
beyond many tens of AU from a solar-type star. 

Any formation theory also needs to explain the wide 
range of physical and orbital properties of giant planets 
around other stars (see the chapter by Cumming). About 
6% of Sun-like stars have giant planets orbiting within 
4 AU. Jupiter-mass planets in orbit within several AU of 
Sun-like stars are more common than planets of several 
Jupiter's masses, and substellar companions more massive 
than ~ 10 Mj are rare. Neptune-mass planets may be 
rather common compared to planets a few times the mass of 
Neptune. These observational results are broadly predicted 
by core accretion models, provided that gas in protoplan- 
etary disks survives for at least a few million years. Disk 
instability models, instead, would predict an abundance of 
massive objects (> 10 Mj) relative to Jupiter-mass planets, 
since formation begins at very early times in massive disks. 
GIs, though, provide an appealing prospect for the forma- 
tion of ~ 10 Mj objects observed at distances > 100 AU 
from their stars, where core formation timescales are prob- 
ably too long. 

Exoplanet data show a strong correlation between stellar 



metallicity and frequency of giant planet detections (e.g., 
Udry and Santos, 2007): ~ 25% of stars whose metallicity 
is twice that of the Sun host a gas giant (within a few AU). 
This percentage reduces to only ~ 5% for solar metallicity 
stars. Disk instability models do not seem to predict this 
trend. In fact, there is some evidence that GIs are stronger 
in lower metallicity disks (Cai et al., 2006). Core accretion 
models do predict this trend (Kornet et al., 2005; Mordasini 
et al., 2009), provided that the density of solids in a disk is 
proportional to the metallicity of the central star. 

Planet searches around M dwarf stars (M* < 0.5 M ) 
seem to suggest that the occurrence of giant planets is less 
likely than it is around more massive stars {Johnson et al., 
2007; Eggenberger and Udry, 2010). Although these sur- 
veys are still incomplete, there is mounting observational 
evidence that planets of mass M < 0.1 Mj orbiting M 
dwarfs may be common (Forveille et al., 2009). If con- 
firmed by more complete statistics, a correlation between 
stellar masses and the masses of hosted planets may pro- 
vide valuable information on the process of planet forma- 
tion. According to core accretion models, gas giants orbit- 
ing low-mass stars should be rare (Kennedy and Kenyon, 
2008), whereas Neptune-mass planets should be common 
(Laughlin et al., 2004). According to disk instability mod- 
els, if conditions are appropriate for disk fragmentation, the 
occurrence rate of giant planets should not depend strongly 
on the stellar mass (Boss, 2006). 

Interpretations of observational data based on either of 
these mechanisms remain difficult and predictions are still 
rather qualitative. There is a number of unresolved issues 
that prevent us from making more quantitative predictions. 
Such issues are generally related to poor constraints on the 
initial conditions used in the formation models. Some of 
these open questions are outlined below. 

What factors determined the final masses of Jupiter and 
other gas giants? The mass of Saturn may have been influ- 
enced by the presence of Jupiter and similar conditions may 
have operated in some extrasolar planetary systems with 
multiple gas giants. However, a more general answer to 
this question will require a much deeper understanding of 
the physical conditions that existed in the solar nebula and 
other protoplanetary disks. 

Are there peaks in the mass distribution function of mas- 
sive planets? During the final phases of growth, giant plan- 
ets are in a symbiotic relationship with their parent disks. 
Therefore, we need better constraints on masses, tempera- 
tures, and lifetimes of protoplanetary disks to make predic- 
tions on the relative abundances of giant planet masses. 

What is the occurrence rate of gas giants and Neptune- 
mass planets as a function of the stellar mass? What is the 
percentage of red dwarfs that host massive planets? Cur- 
rently, a handful of cases are known, as the two Jupiter- 
mass planets orbiting Gliese 876 (http : / / exoplanet . 
eu/), but conclusions cannot be drawn from such small- 
number statistics. 

Is there a numerous population of gas giants orbiting at 
large separations from their host stars, as GIs models seem 
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to predict? If so, how do their orbital properties and enve- 
lope composition differ from those of other known gas gi- 
ants revolving on smaller orbits? What is the likelihood that 
they formed closer to the star (perhaps via core accretion) 
and were later scattered or migrated to larger distances? 

4.2. Impact of Orbital Migration 

Planetary cores of ~ 10 M e that undergo type I migra- 
tion may be subject to rapid orbital decay. One of the main 
advantages of disk instability models is that they would by- 
pass this phase of planetary growth altogether and avoid this 
potential threat. Planets formed via GIs, instead, may be 
subject to prolonged periods of type II migration. But, how 
relevant is planetary migration during planet formation? 

Although a complete answer is not presently known, 
there is evidence that substantial migration may occur dur- 
ing the formation process. A significant fraction of multi- 
planet extrasolar systems are close to or in a mean-motion 
resonance (Udry and Santos, 2007). Since it is unlikely 
that such orbital configurations can be produced by in situ 
formation models, these systems may lend support to con- 
vergent migration hypotheses. Recent studies aimed at ex- 
plaining the structure of the outer solar system (e.g., Mor- 
bidelli et al., 2009) suggest that Jupiter and Saturn may have 
undergone convergent migration, as a result of which Saturn 
became temporarily captured into a mean motion resonance 
with Jupiter. 

The presence of Neptune-mass planets and gas giants or- 
biting stars within a few tenths of AU may also constitute 
evidence of orbital migration, although alternative interpre- 
tations are possible. Planet-planet scattering and Kozai cy- 
cles (see the chapter by Lubow and Ida) can, in principle, 
drive gas giants very close to their host stars. However, 
while these last two mechanisms could produce a large mis- 
alignment between the stellar spin axis and the planet's orbit 
axis, disk-induced orbital migration would preserve spin- 
orbit alignment. Some of the transiting planets allow for 
the measurement of the sky-projected angle between stellar 
spins and orbit axes (see the chapter by Winn). Currently, 
about two-thirds of these measurements are consistent with 
a close spin-orbit alignment (Simpson et al., 2010), suggest- 
ing that migration by tidal interaction with the parent disk 
drove these planets close to their stars. 

The standard theory of orbital migration by disk-planet 
interactions seems able to reproduce the main features of 
the orbital period distribution of exoplanets (e.g.,Armitage, 
2007), which suggest that considerable migration may have 
occurred. Yet, the question of whether or not substantial 
migration is the norm can only be assessed once observa- 
tions of long-period planets become available. A more com- 
plete distribution function of orbital periods can also pro- 
vide quantitative indications on whether current estimates 
of the migration rates of Neptune-mass planets are too rapid 
and, if so, to what extent they are. 



5. FUTURE PROSPECTS 

To achieve a more comprehensive understanding of the 
formation of planetary systems, at least some of the ques- 
tions highlighted in Section 4 need be addressed. We en- 
vision that answers to those questions can indeed be found 
in the near future, as further progress is made along three 
main avenues: data collection from the solar system; obser- 
vations of star-forming regions, protoplanetary disks, and 
exoplanets; theoretical studies of star and planet formation. 

Interplanetary spacecrafts (e.g., Voyager, Galileo, and 
Cassini) have provided a wealth of information about the 
origins and evolution of the outer solar system. The Juno 
mission may finally reveal the mass of Jupiter's core by di- 
rect measurements of the planet's gravity field. Although it 
seems implausible that such information alone will provide 
definitive evidence for one or the other formation scenario, 
it is essential for refining interior models of the gas giant 
nearest to us. 

To better constrain current estimates of sizes, masses, 
dust contents, and lifetimes of protoplanetary disks, it is im- 
portant to know the detailed physical conditions that exist at 
various stages of a disk's evolution, starting from the latter 
phases of stellar formation. The observation of star-forming 
regions and disks surrounding young stars plays a pivotal 
role in establishing more comprehensive and accurate mod- 
els of disk evolution. These models can be then applied 
in the context of increasingly sophisticated calculations of 
planet formation that include disk dispersal. 

Space observatories, like Hubble and Spitzer space 
telescopes, have provided valuable data that has signifi- 
cantly improved our knowledge of star-forming regions and 
planet-forming environments. Over the next decade, a lot 
more is expected to be learned about protoplanetary disks 
with both existing and new-generation instruments, such as 
the Atacama Large Millimeter Array (Tarenghi, 2008), the 
Stratospheric Observatory for Infrared Astronomy airborne 
observatory (Gehrz et al., 2009), and the James Webb Space 
Telescope (http://www.jwst.nasa.gov/). These 
facilities, among others, will improve in completeness the 
mass distribution function of Neptune-mass and larger plan- 
ets. They will likely improve our knowledge on the rela- 
tionships between stellar metallicity and dust content of 
young disks. They should also allow direct imaging of 
high planetary-mass objects and low-mass brown dwarfs, 
orbiting at several tens of AU from stars, and provide pre- 
cious data on the abundances of heavy elements and on the 
chemical compositions of their envelopes. 

Ground-based interferometers (e.g., Very Large Tele- 
scope and Keck telescope) and astrometry space observa- 
tories, such as Gaia (http://sci.esa.int/GAIA) 
and the Space Interferometry Mission (http : //sim. 
jpl.nasa.gov/), will extend the orbital period dis- 
tribution function of detected exoplanets towards longer- 
period orbits. Correlations between planet detection fre- 
quency and stellar metallicity and between stellar mass 
and planet mass are also expected to become more accu- 
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rate and statistically significant. If a significant fraction 
of M dwarf stars possess gas giants, observational data 
should allow us to compare their orbital and physical prop- 
erties with those of planets orbiting more massive stars. If 
a large population of high-mass planets is found at large 
separations, their atmospheric composition can be com- 
pared to those of transiting planets, whose number should 
also grow significantly thanks to space missions such as 
COROT (http://smsc.cnes.fr/COROT/), Kepler 
(http: //kepler. nasa.gov/), and TESS (http: // 
space.mit.edu/TESS). Within the next five to ten 
years, we should be able to test the bimodal distribution 
hypothesis for gas giant formation based on statistical anal- 
ysis of orbital properties, mass distribution functions, and 
envelope compositions. 

The trend of exoplanet detections and characterization 
witnessed over the past decade will carry forward to the 
next decade. Theoretical studies need to tackle still un- 
solved problems on the timescales and physical conditions 
required for gas giant formation. Core accretion models 
must be extended to include proper feedback on the pro- 
toplanetary disk and tidally-induced orbital migration in a 
self-consistent manner. Disk instability models need to be 
refined in order to make quantitative predictions on surviv- 
ability, final orbital radii, and final masses of planetary- 
mass objects formed after disk fragmentation. As data of 
better and better quality become available, it may be chal- 
lenging, from a theoretical standpoint, to combine all pieces 
in a coherent framework. Some of the pieces now regarded 
as important may be discarded and replaced by others. Yet, 
the prospects are in place to turn a blurred picture into one 
with much higher definition within the next ten years. 
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